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ABSTRACT 


A  method  for  passive  ranging  based  on  the  principle  of  triangulation  is  considered 
In  the  basic  triangulation  scheme  that  is  a  single  baseline  model  the  precision  in  the  bearing 
readings  can  be  related  to  the  precision  in  the  range  estimation  For  some  target  orienta¬ 
tions  the  precision  in  the  triangulated  target  range  is  completely  lost  This  phenomenon  is 
known  as  "geometric  dilution."  A  proposed  orthogonal  dual  baseline  scheme  eliminates 
the  "geometric  dilution"  effect.  The  performance  of  each  of  the  two  orthogonal  baselines 
depends  on  target  orientation.  For  specific  target  orientations  the  triangulation  range  mea¬ 
surements  for  the  two  baselines  are  equivalent.  The  dual  baseline  scheme  would  require 
"smart  electronics"  which  would  switch  between  baselines  at  crossover  points  in  the  range 
estimation  precision  It  is  shown  that  the  crossover  points  depend  primarily  on  the  ratio  of 
the  two  baselines.  A  general  expression  for  the  maximum  triangulation  range  consistent 
with  limitations  in  minimum  tolerance  precision  in  range  estimation  is  derived.  The  depen¬ 
dency  between  maximum  range  and  target  orientation  are  presented  in  polar  form.  Limita¬ 
tions  in  the  dual  baseline  model  due  to  physical  limitation  created  by  the  optical  horizon 
are  also  considered. 
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I .  INTRODUCTION 


A.  REVIEW  OF  PRIOR  WORK 

Recent  technological  developments  can  turn  an  active  radar  system  into  a  target  for 
anti-radiation  missiles  (ARM).  In  essence,  a  system  designed  to  provide  security  against 
missile  attacks  can  become  a  security  hazard.  This  thesis  will  present  a  methodology  of 
passive  ranging  based  on  triangulation  and  the  use  of  IR  sensors  which  could  be  imple¬ 
mented  on  various  kinds  of  military  platforms.  IR  sensors  provide  for  passive  surveillance 
and  target  identification  and  tracking,  but  targeting  and  fire  control  require  a  minimum 
range  to  target 

Passive  ranging  is  the  Holy  Grail  of  the  defense  community.  The  passive  interroga¬ 
tion  of  a  scene  greatly  restricts  the  domain  of  applicable  countermeasures.  Therefore,  with 
passive  ranging  the  security  of  the  intelligence  operation  is  less  likely  to  be  compromised. 

Various  passive  ranging  schemes  based  on  the  radiating  characteristics  of  a  target 
have  been  proposed  These  can  be  classified  under  a  variety  of  descriptors  such  as  time 
difference  of  arrival  (TDOA),  frequency  difference  of  arrival  (FDOA),  and  angle  differ¬ 
ence  location  (ADL)  [Ref  1],  [Ref  2],  The  first  two  schemes,  TDOA  and  FDOA  require 
that  RF  modulation  be  present  on  the  target's  electromagnetic  spectrum  [Ref  3]  Since 
the  thrust  of  this  thesis  is  to  focus  on  a  potential  ADL  technique  triangulation,  which  re¬ 
quires  only  a  target's  unmodulated  IR  signature,  no  further  discussion  on  TDOA  or  FDOA 


methods  will  be  included.  Several  methods  for  using  atmospheric  attenuation  to  predict 
target  range  passively  are  covered  in  the  next  paragraph 

In  the  early  1960's  several  patents  were  approved  which  related  IR  signal's  atmos¬ 
pheric  attenuation  to  range.  The  first  of  these  was  based  on  the  atmospheric  absorption  of 
radiation  from  targets  which  have  strong  water  vapor  or  carbon  dioxide  bands  [Ref  4] 
The  scheme  includes  signal  chopping,  optical  beam  splitting,  and  controlled  absorption  of 
the  CO,  band  in  a  reference  path  A  nulling  circuit,  which  uses  a  partially  transmitting 
wedge  driven  by  a  servomotor,  balances  the  two  signals.  The  error  signal  is  then  directly 
related  to  the  amount  of  energy  in  the  CO^  band  of  the  original  IR  signal  received  This  er¬ 
ror  signal  can  be  related  to  the  length  of  the  attenuation  path,  i  .e.,  range.  A  second  patent 
[Ref  5]  is  again  based  on  the  principle  of  preferential  atmospheric  absorption  of  IR  radia¬ 
tion  for  some  bands.  However,  in  this  case  the  authors  proposes  that  if  a  target’s  total  and 
filtered  radiation  intensities  are  known,  and  if  the  effective  radiation  coefficients  for  total 
and  filtered  bands  could  be  ascertained,  then  ratios  of  measured  intensities  can  be  related 
through  analysis  to  the  unknown  range  The  inventors  of  this  patent  also  use  a  chopper 
and  comparison  circuitry.  Recent  NFS  theses  have  examined  the  passive  ranging  problem 
using  similar  concepts  [Ref  6],  [Ref  7],  [Ref  8]. 

The  principle  of  the  parallax  effect,  as  applied  in  the  optical  range  finder  [Ref  9],  has 
been  modified  to  process  IR  signals.  In  this  scheme,  introduced  as  a  patent  [Ref  10],  the 
measurement  of  a  parallax  angle  between  two  imaging  sensors  and  a  known  baseline  will 
permit  calculation  of  the  range  to  a  target.  Another  60's  patent  [Ref  1 1]  suggests  a  model 
for  calculation  closure  time,  the  approximate  time  it  will  take  one  body  to  collide  with 
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another.  This  parameter  depends  on  range.  The  inventors  show  that  the  closure  time  can 
be  calculated  from  the  time  rate  of  change  of  the  target's  passively  detected  radiant  energN' 
They  provide  a  circuit  for  this  application. 

With  the  development  of  IRST  passive  sensors  [Ref  12],  renewed  interest  in  .ADL 
methods  has  occurred.  In  particular  these  devices  are  available  on  a  rotating  (30  rpm)  plat¬ 
form  and  can  predict  bearing  azimuth  data  on  low  flying  targets  A  platform  for  an  IRST 
device  is  shown  in  Figure  1  The  bearing  resolution  of  such  devices  is  classified;  however 
an  interesting  mathematical  model  for  predicting  resolution  for  targets  in  the  presence  of 
noise  has  been  proposed  [Ref  13]. 


Figure  1 :  ERST  sensor  platform  (After  Ref  12) 


j 


Most  of  the  ADL  methods  are  based  on  using  two  sensors,  i  e  ,  a  single  baseline  ap¬ 
proach  (SBL)  [Ref  14],  [Ref  15],  In  the  late  70's  a  multiple  sensor  scheme  based  on  a  te- 
chique  known  as  circulation  [Ref  16],  which  employs  an  interactive  least  squares 
algorithm,  has  been  investigated  Although  the  method  eliminates  the  problem  of  "geomet¬ 
ric  dilution"  the  scheme  can  exhibit  convergence  problems.  More  recently  a  preliminary 
study  on  an  analytic  dual  ba.seline  method  has  been  introduced  [Ref  17]  This  description 
is  based  on  an  approximate  model  for  estimating  the  range  sensitivity  Further  develop¬ 
ments  on  this  dual  baseline  model  will  be  presented  in  this  thesis 

B.  THESIS  OVERVIEW 

Chapter  I  covers  the  literature  review  and  the  thesis  overview.  In  Chapter  II  the  sin¬ 
gle  baseline  analysis  is  developed,  and  the  lack  of  precision  in  a  triangulated  range  is  re¬ 
lated  to  the  lack  <  f  precision  in  the  bearing  measurements.  It  is  shown  that  the  precision  of 
the  triangulation  method  for  targets  alligned  with  the  single  baseline  is  completely  lost 
This  problem  is  known  as  "geometric  dilution  "  Chapter  III  explores  the  applicability  of  a 
dual  baseline  scheme  in  order  to  circumvent  this  problem.  In  Chapter  IV  the  dual  baseline 
model  is  applied  to  characterize  a  maximum  range  for  triangulation  consistent  with  a  mini- 
am  precision  in  the  measurement.  The  development  in  this  chapter  leads  to  a  physical 
polar  representation  for  in  terms  of  target  orientation  Two  paradigms,  quadrature 
addition  and  the  more  conservative  sum  of  absolutes,  are  considered  In  Chapter  V  the 
limitations  in  range  estimation  due  to  earth's  curvature  are  evaluated  and  included  with  the 
dual  baseline  model  Some  interesting  qualifications  for  predicting  the  characteristics  of 
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the  polar  plots  studied  in  Chapter  IV  are  provided  in  Chapter  VI  Chapter  VII  provides 
the  conclusions 


Appendix  A  is  the  list  of  symbols  provided  for  convenience  in  alphabetical  order  In 
order  not  to  complicate  the  main  text  some  of  the  more  mathematical  derivations  have 
been  relegated  to  the  Appendices  B,  C,  and  D  In  Appendix  B  details  on  approximations 
made  on  the  single  baseline  model  are  presented.  In  Appendix  C,  an  exact  mathematical 
expression  or  the  single  baseline  triangulated  range  in  terms  of  IRST  sensors  bearing  is 
derived  In  Appendix  D  two  paradigms  for  error  estimation  are  analyzed.  It  is  shown  that 
the  sum  of  absolutes  estimation  is  more  conservative  than  estimation  via  quadrature  addi¬ 
tion  Appendix  E  is  a  collection  of  the  programs  used  in  generating  the  plots  seen  in  the 
thesis 
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11.  BACKGROUND  ANALYSIS 

A.  THE  TRIANGULATION  CONCEPT 

The  eye-brain  system  undoubtedly  works  in  a  fairly  sophisticated  manner  in  order  to 
extract  information  from  parallax.  The  geometric  considerations  of  the  corresponding  pas¬ 
sive  sensor  problem  are  fairly  simple.  In  principle,  the  range  prediction  can  be  made  with 
1 00%  accuracy;  in  practice  the  limitations  in  the  accuracy  of  the  measurements  in  orienta¬ 
tion  bearings  introduce  an  error  in  the  range  prediction.  The  sensitivity  for  the  range  pre¬ 
diction  to  variation  in  the  bearing  accuracy  is  dependent  on  target  orientation,  the  distance 
between  the  points  of  observation  or  baseline,  and  the  range  to  the  target  Currently  ERST 
bearings  do  not  include  elevation  data.  For  targets  which  appear  coming  over  the  horizon 
this  should  not  represent  a  significant  deficiency  in  the  model.  Therefore,  the  model  to  be 
presented  predicts  horizontal  range,  not  the  slant  range. 

B.  SINGLE  BASELINE  MODEL 

A  single  baseline  model,  involving  two  sensors,  located  at  the  ends  of  a  baseline,  can 
be  analyzed  for  conditions  sufficient  foi  triangulation.  Starting  from  Figure  2  the  points  at 
which  bearing  measurements  are  made,  A  and  C,  are  located  at  a  distance  R,  and  IL  re¬ 
spectively  from  the  target  at  T.  Although  no  measurements  are  made  at  point  B  it  serves 
as  a  convenient  symmetrical  reference  between  observation  points  A  and  C.  The  range  R 
is  defined  as  the  distance  from  point  B  to  the  target.  The  characteristics  of  the  single 
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baseline  model  in  terms  of  symmetric  variables  R  and  0  greatly  facilitates  the  dual  baseline 
analysis  presented  later.  The  distance  between  the  points  A  and  C  is  the  baseline  (BL). 
The  initial  objective  of  this  analysis  will  be  to  relate  the  range  R  to  the  measurables  in  the 
problem,  i.e.,  baseline  (BL)  and  orientation  (0,  and  The  next  step  in  the  procedure  will 
be  to  perform  a  sensitivity  analysis  with  respect  to  variation  in  the  orientation  variables,  0, 
and  03. 

The  significant  assumptions  of  the  model  are; 

( 1 )  There  is  no  error  in  the  baseline  measurement; 

(2)  The  tolerances  in  the  bearing  measurements  are  known,  and  will  be  denoted 
A0,  and  A02 ; 

(3)  The  discrepancy  between  horizontal  path  range  and  slant-path  range  is  not 
significant. 

1.  A  single  baseline  analysis 

From  the  Figure  2,  the  apex  angles  satisfy  the  following  relation 

03  =  03-0,  =  0,  +  0,  (1) 

Another  relation  between  the  apex  angles  and  0,,  0,  is  obtained  from  the  Law  of  Sines 

sin  04  =  ^  sin  02  (2a) 

sin05  =  ^sin0|  (2b) 
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Through  expansion  of  ( 1 )  it  follows  that 


sin(62  -6i)  =  sin64Cos65  +  cos64sin0.<; 


(3a) 


and 


sin  (62-61)  =  ^ 


sin02  Jl  -^^sinOi  j  +sinei  -  ^^sinS- 


(3b) 


Solving  (3)  for  the  range  R,  it  follows  that 


sin  02^ 

fl-^^sin0ij  +sin0)^ 

(l-(fsin02]‘ 

sin(02  -0)) 

(4) 


An  exact  expression  for  R/BL(0,,  0^)  can  be  obtained  using  symbolic  mathematical  proc¬ 
essing.  See  Appendix  C.  The  immediate  objective  is  to  obtain  an  expression  for  the  range 
sensitivity  in  terms  of  the  bearing  angles  0,  and  0^ ,  which  will  be  in  the  form 


^  -/i(01,02)A01  +/2(01,02)A02 


(5) 


In  order  to  not  distract  the  reader,  the  mathematical  details  in  obtaining  the  range  sensitiv¬ 
ity  appear  in  Appendix  B,  and  the  final  expression  for  it  is 
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where  it  is  understood  that  AO,,  and  A6,  are  assumed  to  be  zero  mean  random  variables 
A  statistical  model  for  eliminating  "systematic  error"  has  been  proposed  using  a  scheme 
based  on  a  method  known  as  circulation  [Ref  16],  Furthermore,  the  ratio  "a"  of  rms  val¬ 
ues  for  A0t  and  A0,  can  be  used  and  the  final  expression  for  the  normalized  variation  in  R 
can  be  written  as  follows 


l-||jsin9i  j  {cos02-sin 02cot(0j-ei))-«ui0|  /l-(§ 511182  j  col(03-0|H 

_ 

iinSj  -MuiSi j 


•  X2  •  O 


where  %]  and  Xj  are  independent  random  variables  having  zero  mean  and  unity  standard 
deviation,  and 


(8) 
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It  is  noteworthy  that  (7)  depends  on  R,  0,,  and  6,  instead  of  the  symmetrically  defined 
variables  R  and  6.  Additional  relations  are  needed  at  this  point  in  order  to  express  (7)  and 
(8)  in  terms  of  the  symmetrically  defined  variables  R  and  6.  From  the  Law  of  Cosines, 
applied  to  Figure  2,  it  follows  that 


"i  )'  =  (i)  +(0 


From  the  Law  of  Sines,  applied  also  to  Figure  2,  it  follows  that 


■  n  •  a(«) 

Sin  02  =  sinOjT^ 


(10a) 


sin6i  =sin6“rT 


(10b) 


and,  again  from  the  Law  of  Sines  applied  to  the  same  figure 


sin(02  -0i)  =sin(04  +05)  =  |^sin0i  =  |^sin02 

K\  «2 
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Finally,  along  with  the  trigonometric  rule 


COtX  =  ' 


substitution  of  (8),  (9),  and  (10)  into  (7),  which  is  the  normalized  sensitivity  of  the  range, 
limits  the  dependence  of  variables  to  R  and  6  only.  Although  this  expression  is  not  shown 
explicitly  in  this  form,  it  is  readily  evaluated  in  its  implicit  form  via  the  substitution 
procedure. 

2.  Approximate  single  baseline  analysis  (R»BL) 

To  appreciate  the  factors  that  primarily  limit  the  accuracy  in  the  predicted 
value  for  R  it  can  be  instructive  to  consider  a  simpler  formulation  of  the  problem.  Consis¬ 
tent  with  the  conditions  of  assumption  3  the  ratio  of  the  range  to  the  baseline  can  be  con¬ 
sidered  much  bigger  than  one  (R/BL»1).  After  application  of  this  assumption  to  (3)  the 
second  order  terms  in  BL/2R  can  be  neglected,  leading  from  (4)  to 


p  __  sin j  ^ 

^-2]^  sin{t»2-e,)  J  ~  2 


and  the  normalized  variation  in  range  would  be 


Sin  02-^0101 


■cot(e,-e,)](A9,)+(::;S^-coi(e:-e,))(ae,) 


Equation  (6)  representing  the  normalized  variation  in  R  can  then  be  expressed 


Ag 

K  AH 


sinH:-KinB 


+  cot  (02  -  e , ) )  •  Xi  +  -  cot  (6:  -  e , ) )  •  X:  • 


(15) 


where  again  x,  and  X:  are  independent  random  variables  having  zero  mean  and  unity  stan¬ 
dard  deviation,  and  "a"  is  defined  by  (8).  Equation  (15)  can  re-expressed  in  terms  of  R 
and  0  after  taking  into  account  (9)  through  (12).  Again  the  details  appear  in  Appendix  B 
3.  Approximate  single  baseline  analysis  (R»BL,  R=R,=  R^) 

From  the  Law  of  Sines ,  Equation  (8),  it  follows  that 

8L 

^  1  “  ,/>  ,  sin  01  / 1  z  \ 

sin(e2-«i)  (16a) 


and 


^2=-ifti;ysin02  (16b) 

Working  with  (12),  and  (13)  the  reasonable  assumption  R=R,  =R2  can  be  made  and  (14)  is 
easily  converted  into  the  approximate  form 

^  =  ^^cot0, +cot(02-e,)j(Ae,)-i-(^|cote2-cot(e2-0,)j(Ae:)  (i7) 
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after  noting  that  R  =  R\  =R2  implies 


sui6|  =  sin 02  HsinO 


(18) 


Furthermore  from  (16) 


sin (02  -0i)  =  ^sin6 


(19) 


and  substitution  of  (18),  and  (19)  into  (17)  leads  to 


R 


l-(^sin8) 

f 

1  ,  1 

(Ae,)+ 

1  ^l-sin^0  ^ 

2  sme 

fsme 

2  sm0 

V 

J 

/ 

(Ae.) 


(20) 


Using  the  ratio  "a"  of  A02  and  A0,  applied  to  (19),  it  follows  that 


(  ,  , 

(  ,  1 

l-(^sin9) 

■ 

{ - 

1 V'-**"'*  1 

•Xi  + 

(  ^l-sin^0  J 

R  A9  J 

2  sin0 

&su,9 

2  sin0 

^su)9 

J 

/ 

(21) 


This  expression  represents  the  normalized  sensitivity  in  terms  of  the  correct  symmetrically 
defined  variables,  R  and  0. 

4.  Comparison  of  single  baseline  models 

In  order  to  continue  the  triangulation  analysis  a  comparison  between  the  three 
single  baseline  models.  Equations  (7),  (15),  and  (21),  has  to  be  made  Due  to  very  com¬ 
plicated  expressions  that  stand  for  the  three  models  a  reasonable  comparison  of  them  can 
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complicated  expressions  that  stand  for  the  three  models  a  reasonable  comparison  of  them 
can  be  made  through  graphical  representation.  In  order  to  compare  the  three  approximate 
forms  for  the  range  sensitivity,  the  independence  of  the  random  variables  Xi  and  X-  should 
be  considered.  Since  all  three  approximate  expressions  for  the  range  sensitivity  are  of  the 
form 


TX?  =  Uerml)  •  x,  +  (term!)  •  X: 


(22) 


the  natural  method  of  comparison  is  to  take 

I I  ~  +  (termlf  (23) 

or  equivalently  the  coefficients  of  the  sensitivity  add  in  quadrature.  This,  as  suggested  in 
Appendix  D,  may  not  be  the  most  conservative  estimate  for  the  sensitivity  However  for 
purposes  of  evaluating  the  accuracy  of  the  approximate  model  (21)  relative  to  (7)  and  (15) 
it  is  adequate.  The  three  approximation  levels  of  the  range  sensitivity  expressed  in  Equa¬ 
tions  (7),  (15),  and  (21)  are  compared  in  the  Figures  3,  4,  and  5  for  three  values  of  the 
range  to  baseline  ratio,  and  for  ratio  of  bearing  errors  0,  1  and  2. 

5.  Observations 

Scrutinization  of  Figures  3,  4,  and  5  indicates  only  small  discrepancies  between 
the  approximate  solutions.  It  is  concluded  that  the  approximations  of  Equation  (21), 
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characteristic  of  all  the  single  baseline  figures  are  worth  noting.  First  the  range  sensitivity 
is  highest  for  angles  close  to  0"  and  180“.  This  is  the  well  known  "geometric  dilution"  [Ref 
14]  that  is  characteristic  of  all  the  single  baseline  methods.  In  terms  of  measurements  a 
high  normalized  sensitivity  would  result  in  a  high  level  of  inaccuracy.  This  provides  the 
motivation  for  exploring  the  dual  baseline  method  which  is  the  focus  of  the  later  chapters 
Because  the  range  sensitivity  normalization  is  done  with  the  respect  to  AS,  the  effect  of  in¬ 
creasing  "a"  is  to  amplify  the  term  proportional  to  AG,.  Therefore,  as  seen  from  the  fig¬ 
ures,  there  is  an  overall  increase  in  the  normalized  range  sensitivity.  Lastly,  the  farther  the 
triangulated  target  is  from  the  baseline  the  more  inaccurate  the  measurement.  This  is  also 
confirmed  by  the  increase  in  the  normalized  range  sensitivity  of  Figures  5a,  5b,  and  5c 
over  the  respective  counterparts  of  Figures  3,  and  4. 
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T 


Figure  2 :  Geometry  for  the  single  baseline  model 
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Figure  3a;  Single  baseline  model  sensitivity  curves  for  R/BL=10  &  a=0 
(via  program  in  Appendix  E.  A) 
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comparison  of  single  baseline  sensitivity  curves 
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*  for  more  accurate  case(equation  7) 
o  for  R»BL  approximation  case  (equation  1 5) 

+  for  R»BL  &  R-R1-R2  approximation  case  (equation  21) 
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Figure  3b  :  Single  baseline  model  sensitivity  curves  for  Rj'BL=10  &  a=l 
(via  program  in  Appendix  E.  A) 
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*  for  more  accurate  case(equation  7) 

o  for  R»BL  approximation  case  (equation  1 5) 

+  for  R»BL  &  R-R1-R2  approximation  case  (equation  21) 
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Figure  3c  Single  baseline  model  sensitivity  curves  for  R/BL=10  &  a=2 
(  via  program  in  Appendix  E.  A) 
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comparison  of  single  baseline  model  sensitivity  curves 
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*  for  more  accurate  case(equation  7) 

0  for  R»BL  approximation  case  (equation  1 5) 

+  for  R»BL  &  R-R1  “R2  approximation  case  (equation  21) 
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Figure  4a  .  Single  baseline  model  sensitivity  curves  for  R/BL=100  &  a=0 
(via  program  in  Appendix  E.  A) 
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Figure  4b  Single  baseline  model  sensitivity  curves  for  ft'BL=  1 00  &  1 

(via  program  in  Appendix  E.  A) 
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*  for  more  accurate  case(equation  7) 
o  for  R»BL  approximation  case  (equation  15) 

+  for  R»BL  &  R-R1-R2  approximation  case  (equation  21) 
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Figure  4c  :  Single  baseline  model  sensitivity  curves  for  R/BL=  1 00  &  a= 
(via  program  in  Appendix  E  A) 
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*  for  more  accurate  case(equation  7) 

0  for  R»BL  approximation  case  (equation  1 5) 

+  for  R»BL  &  R-R1*R2  approximation  case  (equation  21) 
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Figure  5a  Single  baseline  model  sensitivity  curves  for  R/BL=  1 000  &  ci=0 
(via  program  in  Appendix  E.  A) 
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Figure  5b  ;  Single  baseline  model  sensitivity  curves  for  R/BL=1000  &  a- 1 
(via  program  in  Appendix  E.  A) 


24 


Idelta  R  /  (R  •  delta  theta)| 


xIO 

2.5| - r 


cxsmparison  of  single  baseline  model  sensitivity  curves 


*  for  more  accurate  case(equation  7) 
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Figure  5c  :  Single  baseline  model  sensitivity  curves  for  R/BL=1000  &  a=2 
(via  program  in  Appendix  E.  A) 
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m.  DUAL  BASELINE  MODEL 

A.  DUAL  BASELINE  MODEL  ANALYSIS 

Figure  6  shows  a  dual  baseline  model.  The  two  baselines  of  the  model  are  taken  to 
be  orthogonal  and  they  are  referred  to  as  the  long  base  line  (LBL)  and  the  short  baseline 
(SBL)  system  The  features  and  the  relations  of  the  single  baseline  model  analysis  are  ap¬ 
plicable  here  but  the  superscripts  S  and  L  will  be  used  to  distinguish  the  two  systems. 
Furthermore,  some  new  relations  between  the  angles  are  needed,  in  order  to  express  the 
short  baseline  (SBL)  model  in  terms  of  the  same  variables  R  and  6,  as  in  the  long  baseline 
(LBL)  model. 

B.  GENERAL  CONSIDERATIONS 
From  the  Figure  6  it  is  clear  that 

cos  0^  =  cos  (6  -  90’)  =  sin  0  (24) 


similarly 


sin0^  =  -cos0 


Therefore  from  Equations  (9),  (24),  and  (25) 

[m]  =[m]  +(i)+(2lz:)cos0 


(25) 


(26a) 
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(26b) 


sfii)  -(ir)  +(2)  +(m)s>ne 


(26c) 


(26d) 


Similarly,  from  Equations  ( 1 0),  (24),  and  (25)  it  follows  that 


sin  of  =  sin  0^^^ 


f-£- 

sin  62  =  sin  67^ 


(27a,  b) 


sin  0]  = 


-COS0 


V«n  sin0f  =  -cos0^ 


(27c,  d) 


and  from  Equations  (11),  (24),  and  (25) 


sin  ^0^  -  0t  j  =  ^  sin  ^  sin  0^ 


(28a) 


sin  [ef  -  0f  j  =  ^  sin  0|  =  ^  sin  0^ 


(28b) 
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For  computation  of  j 

(S:)~^(z^]  (S)^(m]  01-^01  02^02  (29a,  b,  c,  d) 


and  follow  substitution  (26a,  b)  and  (27a,  b).  For  computation  of 


f  ~  t' 


let 


01— >0‘^  02  — >  02  (30a,  b,  c,  d) 


and  follow  substitutions  (26c,  d)  and  (27c,  d).  In  Appendix  B  expressions  for  j 
and  j  are  shown  for  the  various  levels  of  approximation  introduced  in  the  previ¬ 

ous  single  baseline  description.  As  noted  in  this  chapter  the  approximations  R»BL,  and 
Ri=R,=R  are  quite  good,  and  the  dual  baseline  results  (see  (21))  for  this  case  are  provided 
here  for  later  reference 


( 

■ 

f 

(  _ 

\  \ 

J  ^ l-sm^6  1 

(  /?  ae  J  ~ 

2  sin  6 

> 

Xi 

2  sin  9 

(¥•>"») 

•X2«'-(31a) 


- 

f  h  ,,nV 

( 

fi-(^cose] 

]  Jl-cm^6  ij  {  >>  ■  j 

]  ^i-cos^e  jj 

[  R  A8  J  ~ 

2  (-cose)  ' 

•  Xi  + 

2  (-COS0) 

\ 

•X^«''■(31b) 
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where  X\  '  Xz  ^  X:^  3re  independent  random  variables  having  zero  mean  and  unity 

standard  deviation,  and 


= 


a^  = 


(32a,  b) 


C.  SENSITIVITY  CURVES  FOR  THE  DUAL  BASELINE  SCHEME 

The  expressions^— ^  j  ,  and  j  in  (31)  represent  the  normalized  sensitivity 
for  the  long  and  the  short  baseline  system  respectively,  in  terms  of  the  symmetrical  de¬ 
pendent  variables,  R  and  0.  It  is  obvious  from  (31)  that  the  sensitivities  of  both  long  and 
short  baseline  systems  depend  heavily  on  the  terms  {R/LBL),  {R/SBL\  AO,'-,  AB^^  AG,", 
and  AG/  To  illustrate  the  dependence  of  the  sensitivities  to  the  above  mentioned  terms 
some  examples  will  be  introduced  with  the  use  of  (31)  which,  after  (23),  becomes 


( 

1  y  l-sin^0  ) 

r.L 

2  sin0 

( LBL  ■  rs  ^ 

•  a 

nr*'"® 

. 

V  «  }  ) 

(33a) 


1  Vi-cos^e  y 

2  (-COS0)  I  SSL 


2  \ 


I  » 


(-cos 


e,) 


J_ _ 

2  (-cos  0) 


LV 


(  SSL 
V  K 


(-COS0) 


■  a’ 


(33b) 
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where  as  previously  described  the  terms  are  added  in  quadrature  However  geometrical 
considerations,  described  in  Appendix  D,  suggest  that  a  more  conservative  gauge  for 
range  sensitivity  is  based  on  addition  of  absolutes 

I  =  Werm\)\  +  \iterm2)\  (34) 

under  this  assumption  Equation  (3 1 )  would  be  evaluated  according  to 


( 

+ 

]  ^l-sin^0 

( Ae  j 

2  sin  6 

2  sin  6 

(35a) 


f  AR  1 

-t- 

(  «  A9  j 

2  Mose)  '  (^^(^0,9)) 

2  (^OSB)  (^(^ose)) 

(35b) 


It  can  be  seen  from  Appendix  D  that  for  the  discussion  given,  (35)  for  r/  =  cf=  1  is  nomi¬ 
nally  more  conservative  than  (33).  Furthermore  for  c^=  0  ,  Equation  (33)  and  Equa¬ 

tion  (35)  are  equivalent. 

D.  EXAMPLES 

The  examples  discussed  in  this  chapter  will  be  based  on  the  sensitivity  expression  us¬ 
ing  quadrature  addition.  Equation  (33).  The  concepts  applied  are  equally  applicable  to  the 
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more  conservative  equation  (35).  MATLAB  programs  for  evaluating  (33 

Appendix  E  As  a  first  example  take 

)  are  provided  in 

^  =  600 

(36a,  b) 

=2 

(36c,  d) 

After  substitution  of  the  data  into  (33)  and  plotting  it  the  results  are  shown  on  Figure  7. 

As  a  second  example  take 

^=1000 

Wi  =  5000 

(37a,  b) 

a^  =  1 

a^=l 

(37c,  d) 

The  plots  following  substitution  of  (36),  and  (37)  into  (33)  are  shown  in  Figure  8.  In 
both  Figures  7  and  8  the  performance  of  the  short  baseline  system  will  exceed  that  of  the 
long  baseline  for  angles  close  to  0°  and  1 80”.  For  Figures  7  and  8  those  crossing  points  are 
approximately  +17”  and  +15",  respectively.  To  appreciate  the  utility  of  the  analysis  take  a 
simple  example  (i.e.,  that  one  sensor  on  each  baseline  has  no  bearing  error, while  the  bear¬ 
ing  errors  of  the  remaining  sensors  are  equal).  Then  the  situation  can  be  described  as 
follows: 

A0j=Aef  =  Ae  (38a) 
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ABj  =  ABj  =  0 


(38b) 


thus,  the  long  and  the  short  parts  of  (33)  after  setting  =  0  become 


M_L  _ 

/f  46  J  2  sine 


J]-[^s,n9) 


(39a) 


M_L  'l  I  i/<-co^ 

/?  46  j  2  cose 


^  Ji-{^cose) 


-cose 


(39b) 


which  permit  the  plotting  of  the  sensitivities  versus  B.  This  was  done  for  the  conditions 

^=  1000  and  ^  =  5000  (40a,  b) 

and  the  results  are  represented  on  Figure  9.  Note,  as  expected,  the  LBL  does  better  over 
a  wider  range  of  angles.  Around  0=170”  and  B=10”  there  is  a  crossover  point  for  which 
the  short  baseline  system  would  produce  more  accurate  range  measurements  at  angles 
B<10'’  and  B>170'’  For  illustration  of  the  use  of  the  curves  assume  that  the  bearing  resolu¬ 
tion  is  AB‘'|  =A6^2  =AB  =  0.1  mr  and  the  range  to  the  target  is  20  Km.  Note  according  to 
(39)  this  implies  long  and  short  baselines  of  20  m  and  4  m,  respectively  For  B=90°,  the 
LBL  best  case,  the  uncertainty  in  the  long  baseline  range  measurement  from  (39a)  is 

=/?•  ABtanv|/^  (41) 
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where  AR^  is  the  uncertainty  in  the  long  baseline  range,  and 


\|/^  =atan 

Then,  the  range  measured  by  the  long  baseline  will  be 

RLas=R±^ieas  (43a) 

R^„^=20  Km  ±  2.0  Km  (LBL)  (43b) 

For  6=0“  (or  1 80°),  the  SBL  best  case,  the  uncertainty  in  the  short  baseline  range  measure¬ 
ment  from  equation  (39b)  is 

^ieas  =  ^  •  Ae  •  tan  (44) 

where  AR®  is  the  uncertainty  in  the  short  baseline  range,  and 


Thus,  the  range  measured  by  the  short  baseline  will  be 


(45) 


Rieas  =R±ARLos 


(46a) 


R'me«  20  Km  ±]0  Km  (SBL) 


(46b) 


Although  the  inaccuracy  in  the  measured  value  for  R  may  be  tolerable  in  the  LBL  case,  the 
SBL  range  estimate  is  probably  not  acceptable  Figure  9  shows  what  happens  to  the  nor¬ 
malized  range  sensitivity  versus  6  for  fixed  range  For  the  dual  baseline  scheme  the  meas¬ 
urement  accuracy  is  at  a  minimum  at  the  crossover  points.  For  the  conditions  taken  this 
occurs  at  6=+ 10“  and  6=-i-170“. 

E.  RANGE  LIMITATIONS  FOR  FIXED  TARGET  ORIENTATION 

Rang"!  sensitivity  curves  plotted  versus  range  for  the  principal  directions  6  =  0“  and  6 
=  90“  could  be  applied  to  predict  maximum  bearing  error  and/or  minimum  baselines  in  or¬ 
der  to  accurately  triangulate  targets  at  desired  distance  limits.  Under  the  assumptions  in¬ 
troduced  by  (38)  the  normalized  range  sensitivities  are  given  in  (39),  which  for  the  single 
baseline  model  becomes 


1  J l-sin^(8) 

2  sin  (6) 


1(0)1 


(fsm(e)) 


(47) 


Equation  (47)  is  plotted  with  R/'BL=\00  to  R/BL=\000,  and  for  angle  0=90",  45",and  22  5" 
on  Figure  10.  In  order  to  clarify  the  application  of  Figure  10  to  dual  baseline  triangulation, 
an  example  is  presented  based  on  the  following  conditions 

LBL=20.0  m  SBL=4.0  m  (48a,  b) 

0=90“  A0=O.  1  mr  (48c,  d) 
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and  precision  in  range  <  5%,  that  means 


f  <  0.05  (48e) 

From  (48d),  and  (48e)  the  range  sensitivity  is  evaluated  as 

(fi)<500  (49) 

and  therefore  from  Figure  10  the  range  to  baseline  ratio  is  found  to  be 

K  S  550  (50) 

This  condition  is  applicable  to  both  the  long  and  short  baseline  In  both  cases  the  90"  con¬ 
dition  (48c)  imposes  the  interpretation  that  the  target  appears  along  a  line  which  is  perpen¬ 
dicular  to  the  baseline. 

For  the  LBL  case  it  follows  from  Equation  (48a)  that 

^  550x20 w  =  1  IAtw  (51a) 

and  similarly  for  the  SBL  case 

Rsbl  X  4.0  m  =  2.2  km  (51b) 
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This  demonstrates  that  the  target  must  be  significantly  closer  to  be  triangulated  with  5% 
accuracy  if  the  oncoming  direction  is  perpendicular  to  the  short  baseline  Nonetheless,  as 
previously  mentioned,  the  closest  unresolved  range  will  occur  for  the  previously  noted 
crossover  angle.  To  demonstrate  the  relevance  of  this  issue.  Figure  1 1  shows  the  plot  of 
(47)  for  the  angles  0=90”,  0=8O”,and  0=10”.  The  angle  10°  was  chosen  because  it  repre¬ 
sents  the  approximate  crossover  point  for  the  conditions  (38).  After  referring  to  Figure  1 1 
it  is  seen  that  the  performance  for  the  0=10°  case  is  about  six  times  worse  than  the  0=90" 
case.  Using  the  LBL  baseline  (51a)  the  range  detected  by  the  long  baseline  would  be 
approximately 


R<\.^km  (52) 

which  is,  as  expected,  worse  than  either  of  the  baseline  performances  (51).  It  should  be 
noted  that  for  a  crossover  point  the  long  base  line  and  short  baseline  performances  are  the 
same.  This  can  be  checked  by  using  the  0=80”  curve  with  the  SBL  data. 

F.  COMMENTS  ON  THE  EXAMPLES  AND  ANALYSIS 

The  primary  objective  of  the  preceding  examples  and  discussion  has  been  to  demon¬ 
strate  a  feasible  modeling  scheme  for  dual  baseline  triangulation  The  results  of  this  chap¬ 
ter  demonstrate  several  features  of  the  dual  baseline  scheme.  The  problem  of  "geometrical 
dilution"  discussed  in  Chapter  I  can  be  avoided  by  using  two  orthogonal  baselines  In  gen¬ 
eral  the  larger  baseline  will  have  a  higher  measurement  accuracy.  Because  of  typical  ship 
structure  it  is  assumed  that  a  second  baseline  would  not  be  conveniently  as  long  The 
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scheme  would  require  some  "smart  electronics"  which  would  switch  between  baselines  at 
a  performance  crossover  point.  As  seen  from  the  discussion  in  section  D  the  maximum  tri¬ 
angulation  range,  for  a  fixed  range  sensitivity  ,  can  be  predicted  for  particular  target  ori¬ 
entations,  i.e.,  fixed  0.  In  the  chapter  to  follow  this  concept  is  extended  by  solving  for 
versus  0  under  conditions  of  constant  range  sensitivity. 
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Figure  8:  Dual  baseline  second  example(R/LBL=  1 000,  R/SBL=5000,  a^=l,  a"=l) 
(via  program  in  .Appendix  E.B) 
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Figure  9:  Dual  baseline  third  example(R/LBL=1000.  R/SBL=5000,  ct -£^=0) 
(via  program  in  Appendix  E.  B) 
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Figure  10;  Range  sensitivity  for  fixed  target  onentation  (0=  90°,  45°,  and  22.5°) 
(via  program  in  Appendix  E  C) 
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Figure  1 1 :  Range  sensitivity  for  fixed  target  orientation  (6=  90“.  80°,  and  10“) 
(via  program  in  Appendix  E.  C) 


IV.  RANGE  LIMITATIONS  FOR  FIXED  SENSITIVITY 

A.  CHAPTER  OVERVIEW 

In  this  chapter  a  more  graphical  representation  for  versus  orientation  angle  6 
will  be  developed.  For  target  ranges  less  than  the  triangulation  estimation  will,  by  as¬ 
sumption,  be  within  a  tolerable  precision  limit  dictated  by  a  maximum  sensitivity.  Two 
paradigms  for  estimation  are  considered.  The  first  is  based  on  a  quadrature  addition  of 
terms.  In  this  case  an  analytic  solution  can  be  obtained  for  the  case  a=\,  i.e.,  equal  preci¬ 
sion  bearing  measurements.  In  the  second  case  the  addition  occurs  as  a  sum  of  absolutes. 
Here  after  making  a  reasonable  approximation  R/BL  »1,  a  solution  for  all  "a"  can  be 
obtained. 

B.  ANALYSIS  BASED  ON  THE  SUM  OF  ABSOLUTES 

Define  the  normalized  range  sensitivity  for  long  and  short  baseline  respectively  as 


^_L 
R  AB 


(53a,  b) 


after  noting  for  both  LBL  and  SBL  that 


2  iin(8)  ^sui(e) 


(53c) 
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and  after  substituting,  (35)  simplifies  to 


(54a) 


(54b) 


For  typical  cases  R  »BL  (as  previously  assumed  in  Chapter  II  )  the  first  terms,  propor¬ 
tional  to  (1-0^^),  are  very  small  compared  with  the  second  ones,  proportional  to  (l+o^'*), 
so  it  is  possible  to  neglect  them  without  significant  sacrifice  in  the  accuracy.  Then  (54)  be¬ 


comes 


(^$111(64)) 


(1  +a^)  =  c^ 


(55a) 


(^»ir(9^)) 


(l-l-o^)  =  c^ 


(55b) 


Solving  (55a)  for  it  follows  that 


/  \  2 

[Srj  -sin  (6)-^ 


LBL  ~  t 


(55c,  d) 
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Equation  (55d)  can  be  used  to  create  a  polar  plot  of  the  range  measured  by  the  long  base¬ 
line  normalized  to  the  LBL  Solving  (55b)  for  SBL,  and  after  (25),  it  follows  that 


f<SBl 

SBL 


cos‘(6  ) 


(l+a-'r 


(55e) 


Bsbl 

LBL 


COS*(0) 


(55f) 


Equation  (55f)  gives  the  range  measured  by  the  short  baseline  normalized  with  respect  to 
the  LBL  for  convenient  comparison  with  (54). 

C  ANALYSIS  BASED  ON  QUADRATURE  ADDITION  OF  TERMS 

After  the  reasonable  assumption  that  a  =  1,  (33)  becomes 


f 

fAff  1  '1^  _ 

*  sui(9^) 

V 

T 

2  sin(0^) 

y 

(56a) 


fAR  1  1 

2 

f 

-j 

2  ™(Bh  ' 

T 

2  sin  (9*) 

k 

(56b) 


After  defining 


1 

2  sin(H) 


COt(9)  =:^ 


(57) 
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and 
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After  referring  to  Figure  5  and  solving  (60a)  for  it  is  found  that 


f<LBL 

LBL 


fi 


)— 0  5cos^(6 


+  sin‘(6  ) 


(61a) 


Similarly,  solving  (60b)  for  and  taken  into  account  (24)  it  follows  that 


fisBL 

LBL 


(c^)"cos^(Q)-0  5sin‘(8) 


(61b) 


Equations  (6 1  a,  b)  can  be  plotted  in  polar  form  in  order  to  show  the  target's  range  (nor¬ 
malized  to  LBL)  as  it  would  be  measured  by  a  long  and  short  baseline  respectively 

D.  RESULTS  BASED  ON  SUM  OF  THE  ABSOLUTES 

Following  (55d)  and  (55f),  the  maximum  normalized  range  LBL  and  SBL  have  been 
plotted  in  Figures  12,  and  13  respectively,  versus  target  orientation  6.  T’nis  was  done  for 
conditions  a^=a^=l,  c^=c^=500,  and  LBL=5  x  SBL.  For  example,  using  0=90°  on  Figure 
1 2  indicates  that  the  maximum  triangulation  range  is  250  x  LBL.  It  is  seen  from  Figure  1 3 
that  the  SBL  maximum  occurs  at  0“  and  it  is  50  x  LBL  as  expected.  The  combined  effect 
is  easily  visualized  in  the  cartesian  representation  shown  on  Figure  14.  The  solid  line  is  the 
maximum  triangulation  range  using  the  short  baseline  while  the  small  circle  plotting  sym¬ 
bol  is  the  corresponding  long  baseline  range.  As  mentioned,  the  assumption  inherent  in  this 
scheme  is  that  the  baseline  producing  the  most  precise  estimate  for  range  will  automati¬ 
cally  be  selected  with  the  use  of  "smart  electronics."  Therefore  the  larger  of  the  short  and 
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long  baseline  curves  for  represents  the  ideal  capabilities  for  the  dual  baseline  scheme 
This  is  shown  on  Figure  14  with  a  large  circle  used  as  a  plotting  symbol  The  more  physi¬ 
cal  polar  representation  for  this  solution  is  provided  in  Figure  15.  In  this  figure  can  be  seen 
the  maximum  range  which  is  measured  via  the  dual  baseline  scheme  in  any  direction. 

E.  RESULTS  BASED  ON  QUADRATURE  ADDITION  OF  TERMS 

Figures  16,  17,  18,  and  19  are  the  counterparts  for  Figures  12,  13,  14,  and  15  with 
the  same  data  (i.e.,  a^=a^=\,  c^==c^=500,  and  LBL=5  x  SBL)  but  based  on  {61a),  and 
(61b).  It  should  be  noted  that  these  curves  have  a  that  as  expected  is  approximately 
2'  ■  higher  (i.e.,  less  conservative)  compared  with  the  previous  case. 

F.  OBSERVATIONS 

It  is  possible  to  obtain  approximate  analytic  expressions  for  /  LBL  in  terms  of 
the  target  orientation.  In  the  case  of  the  more  conservative  sum  of  absolutes  this  can  be 
done  for  arbitrary  bearing  precision  ratios  and  (f.  Figure  12,  obtained  using  the  sum  of 
absolutes  paradigm,  shows  that  the  maximum  value  for  the  long  baseline  range  is  approxi¬ 
mately  250  X  LBL,  where  250  is  50%  of  the  long  baseline  sensitivity  This  can  be  com¬ 
pared  with  similar  long  baseline  results  for  the  quadrature  addition  paradigm  shown  on  the 
Figure  16.  Here  the  maximum  value  for  the  long  baseline  range  is  approximately  350  x 
LBL.  where  350  is  the  70%  of  the  long  baseline  sensitivity.  These  observations,  which 
show  that  the  range  scales  with  the  sensitivity,  can  be  generalized  to  convenient  rules  of 
thumb  for  calculating  the  maximum  triangulated  range  for  the  long  baseline  system  The 
details  of  this  case  are  presented  in  Chapter  VI.  This  suggest  that  for  the  condition  c‘ 
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c  c  the  general  appearance  of  the  two  curves,  independent  of  scale,  will  depend  on  the 
ratio  of  the  long  baseline  to  the  short  baseline.  This  issue  is  clarified  also  in  greater  depth 
in  Chapter  VI. 
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normalized  range  detected  by  long  baseline  (sum  of  absolutes) 
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ratio  of  the  bearing  errors  in  long  baseline=1  normalized  sensitivity  of  the  long  baseline  =500 


Figure  12:  Normalized  range  measured  by  the  long  baseline  sensors  for  a^=l,  and  c’^=500 
(based  on  the  sum  of  absolutes  and  obtained  via  program  in  Appendix  E.  D) 
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normalized  range  detected  by  the  short  baseline(sum  of  absolutes) 


ratio  of  bearing  errors  in  short  baselin=1  normalized  sensitivity  of  the  short  baseline=500 


Figure  13:  Normalized  range  measured  by  the  short  baseline  for  a^=l.  c^=500,  and 
SBL/LBL=0.2  (based  on  the  sum  of  absolutes  and  obtained  via  program 
in  Appendix  E.  E) 
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range  detected  by  the  dual  baseline  system  in  any  direction(sum  of  absolutes) 


Figure  14:  Cartesian  plot  of  the  normalized  range  measured  by  the  short  baseline ,  and  the 
long  baseline,  showing  also  the  range  measured  by  the  entire  dual  baseline  sys-- 
stem,  (based  on  the  sum  of  absolutes  and  obtained  via  program  in  Appendix 
E.  F) 
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range  detected  by  the  dual  baseline  system  in  any  direction(sum  of  absolutes) 
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ratio  of  bearing  errors  in  both  baselines=1  normalized  sensitivity  in  both  baselines=500 


Figure  15  Plot  of  the  normalized  range  measured  by  the  dual  baseline  system  in  any  di¬ 
rection.  (based  on  the  sum  of  absolutes  and  obtained  via  program  in  Appen¬ 
dix  E.  G) 
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normalized  range  detected  by  long  baseline(quadrature  addition  of  terms) 


270 

ratio  of  the  bearing  errors  in  long  baseline=1  normalized  sensitivity  of  the  long  baseline=500 


Figure  16  Normalized  range  measured  by  the  long  baseline  sensors  for  and  c‘'=500 
(based  on  the  quadrature  addition  &  obtained  via  program  in  Appendix  E  H) 
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normalized  range  detected  by  the  short  baseline{quadrature  addition  of  terms) 


270 

ratio  of  beanng  errors  in  short  baseline=1  normalized  sensitivity  in  short  baseline=500 


Figure  17:  Normalized  range  measured  by  the  shon  baseline  for  a^=l,  c’^=500,  and 

SBL/LBL=0,2  (based  on  the  quadrature  addition  and  obtained  via  program 
in  Appendix  E.  I) 
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theta  in  degrees 


Figure  18:  Cartesian  plot  of  the  normalized  range  measured  by  the  short  baseline  .  and  the 
long  baseline,  showing  also  the  range  measured  by  the  entire  dual  baseline  sy¬ 
stem  (based  on  the  quadrature  addition  and  obtained  via  procram  in  Appendix 
E  J) 
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range  detected  by  the  dual  baseline  system  in  any  direction(quadrature  addition) 
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ratio  of  bearing  emors  in  both  baselines*1  normalized  sensitivity  in  both  baseiines«500 


Figure  19:  Plot  of  the  normalized  range  measured  by  the  dual  baseline  system  in  any  di¬ 
rection.  (based  on  the  quadrature  addition  and  obtained  via  program  in  Appen¬ 
dix  E.  K) 
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V.  PHYSICAL  LIMITATIONS  OF  THE  TRIANGULATION  METHOD 


A.  INTRODUCTION 

The  polar  plots  of  the  normalized  sensitivity  that  were  presented  in  the  previous  sec¬ 
tions  define  limitations  in  the  range  measurement  via  the  triangulation  model.  In  this  chap¬ 
ter  the  physical  limitation  due  to  the  optical  horizon  is  considered. 

B.  GEOMETRIC  AND  OPTICAL  HORIZON  [REF.  18) 

The  line  at  which  the  earth's  curved  surface  apparently  meets  the  sky  is  called  the  ho¬ 
rizon  It  can  be  observed  very  easily  that  the  distance  to  the  horizon  increases  as  the  height 
of  an  observer's  eye  increases  above  the  earth  surface.  The  dependence  of  sighting  dis¬ 
tance  on  the  height  of  the  observer's  eye  is  by  itself  an  evidence  that  the  surface  of  earth  is 
curved  As  shown  in  Figure  20,  two  types  of  horizons  can  be  distinguished:  the  geometric 
horizon  that  extends  to  the  point  H,  and  the  optical  horizon  that  extends  to  the  point  //'. 
The  geometric  horizon  can  be  defined  as  the  points  where  all  the  straight  lines  from  a 
given  observer  tangentially  graze  the  surface  of  the  earth.  However,  due  to  refraction  an 
observer  can  see  a  longer  distance  than  that  of  the  geometric  horizon.  This  displaced  hori¬ 
zon  is  the  only  one  that  can  be  sensed  by  an  observer  and  it  is  called  the  optical  horizon 
So,  the  refraction  enables  the  observer  to  see  over  and  beyond  the  earth's  surface  at  H. 
The  distance  to  which  an  observer  can  see  depends  on  the  height  of  the  eye  and  the 
amount  of  refraction,  and  uncertainties  in  this  latter  quantity  can  be  troublesome  when 
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high  accuracy  is  needed.  The  uncertainties  arise  from  the  fact  that  the  line  of  sight  to  the 
horizon  lies  too  close  to  the  earth's  surface.  Because  of  differences  in  air  and  land  or  water 
temperature,  abnormal  temperature  inversions  or  temperature  gradients  may  occur  in  the 
lowest  few  meters  of  the  atmosphere.  Thus,  the  amount  that  the  line  of  sight  is  bent  is  a 
variable  and  generally  very  difficult  to  predict,  leading  to  a  high  level  unpredictability  in 
the  optical  horizon  location.  Therefore  the  more  conservative  estimation  for  the  optical 
horizon,  based  on  geometric.il  consideration  alone,  will  be  followed. 

C.  OPTICAL  HORIZON  GEOMETRIC  ANALYSIS 

As  may  be  understood  from  the  Figure  20,  no  part  of  the  earth's  surface  whose  ele¬ 
vation  angle  is  less  than  that  of  the  optical  horizon  can  be  seen  beyond  this  horizon  How¬ 
ever,  an  elevated  object  such  as  an  incoming  missile,  flying  above  the  surface  of  the  sea, 
can  be  seen  beyond  this  horizon.  The  geometry  for  a  such  an  object  is  shown  on  Figure 
2 1 ,  and  the  foregoing  analysis  refers  to  an  average  steady-state  condition  of  the  air.  From 
Figure  2 1 ,  it  is  obvious  that 


d  =  {Re  + 1)  •  sin  6,  +  {Re  +  5)  •  sin  0 j 


(65) 


Again  from  the  same  figure 


COS0, 


R,-h 


(66a) 


(66b) 
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After  (66),  (65)  becomes 


c/ =(/(.  +  /) 


+  (Re  +-i') 


(67a) 


or 


d=  J{t-  +  2Rei)  +  J(s-  +  IRfS) 


(67b) 


In  (67b)  the  terms  r  and  r  are  very  small  compared  to  the  other  terms  wiln  factor  and 
can  be  neglected  without  error.  After  this  assumption  (67b)  becomes 


dB  JlRet  +  JlReS  (68) 

The  geometric  analysis  ignores  the  effects  of  the  refraction  which  can  increase  the  distance 
to  the  optical  horizon.  Therefore  (68)  represents  a  conservative  estimation  for  the  optical 
horizon  The  relation  (68)  gives  an  approximation  of  the  distance  at  which  an  object  at 
height  /  above  the  surface  of  the  earth  can  be  seen  by  an  observer  or  sensor  at  height  i 
above  the  surface  of  the  earth.  Taking  R^=6, 37x10  *  meters.  (68)  can  be  written  as 

=  3569.6(  y?  +  .yr )  meters  (69) 

where  d,  t .  and  s  are  expressed  in  meters. 
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D.  EXAMPLES  INCORPORATING  THE  GEOMETRIC  HORIZON 


The  dual  baseline  triangulation  model  permits  a  maximum  triangulation  range  to  be 
estimated  within  a  maximum  tolerable  error  criterion.  In  the  case  that  the  geometric  hori¬ 
zon  appears  at  a  distance  shoaer  than  the  dual  baseline  prediction  then  the  de-facto  maxi¬ 
mum  triangulation  range  is  the  geometric  horizon.  The  combined  limitations  of  the  dual 
baseline  triangulation  system  composed  of  the  inherent  and  physical  limitations  can  be 
modeled  and  represented  with  computer  plots.  Such  plots  are  presented  in  Figures  (22a, 
b)  and  (23a,  b)  and  it  can  be  seen  that  the  limitations  in  the  behavior  of  the  system  depends 
either  on  the  distance  from  the  baseline  to  the  geometric  horizon  or  the  parameters  of  the 
system  itself  The  data  used  to  generate  these  plots  are  conveniently  provided  in  Table  1 . 
It  is  noted  that  the  specification  </=</=!  is  consistent  with  the  criterion  that  all  bearing 
sensors  have  the  same  precision  in  bearing  measurement.  On  the  other  hand  the  specifica¬ 
tion  6-^=  c‘^=c  is  consistent  with  the  condition  that  the  range  measurement  should  be 
worked  with  the  same  precision  limit  regardless  of  the  baseline  used.  Note  that  Figures 
22a  and  23a  are  based  on  the  sum  of  absolutes,  while  the  Figures  22b  and  73b  are  based 
on  the  quadrature  addition  of  terms  On  all  figures,  a* ** ***'=a^=l,  and  c‘^=6'^=500. 

TABLE  1  DATA  FOR  THE  FIGURES  (22,  23) 


t(m) 


s  (m) 


GH  (m) 


LBL  (m) 


SBL  (m) 


Figures  22a,  23a 
Figures  22b,  23b 


20 

3 


10  22,480 

3  12,365 


100 

20 


*  On  all  figures  =0^=1  ,  and  c^=c‘'’=500 

**  Figures  (22a),  and  (22b)  are  based  on  the  sum  of  absolutes 

***  Figures  (23a),  and  (23b)  are  based  on  the  quadrature  addition  of  terms 
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Figure  20  Horizon  geometry  and  efTects  of  refraction  (After  Ref  18) 
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Figure  2 1  Geometrical  horizon  geometry  in  case  that  both  sensor  and  target  are  above 
the  surface  of  the  earth 
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detected  range  is  limited  by  the  geometric  horizon  (sum  of  absolutes) 


90 

2.5e+004 


270 


height  of  the  target  above  surtace^lOm  height  of  the  sensor  above  surf  ace*  10m 


Figure  Z2a.  The  maximum  measured  range  is  iimiied  due  co  geomerric  horizon  in  22A80 
m  instead  of  25. 000m  that  can  be  detected  by  the  system,  (based  on  the  sum 
of  absolutes  and  obtained  via  program  in  .Appendix  E.  L) 


d@tected  range  is  limited  by  the  system  (sum  of  absolutes) 

90 

1.4e+004 


range  detected  by  the  system 


geometric  horizon 


180 


long  baseline  length-20m 
short  baseline  length*4m  240 

height  of  the  target  above  surtace-3m 


270 


height  of  the  sensor  above  surface* 


Figure  Z2b  The  maximum  deteaed  range  is  limited  to  5.000m  due  co  rhe  method,  instead 
the  I2.j65m  which  is  the  geometnc  horizon  (based  on  tne  sum  of  absolutes 
and  obtained  via  program  in  .Appendix  E.L  ) 


'3m 
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detected  range  is  limited  by  the  geometric  horizon  (quadrature  addition) 


®°2.5e+004 


Figure  23  a  The  maximum  measured  range  is  limited  due  to  geometric  horizon  in  22,480 
m  instead  of  3 5,000m  that  can  be  detected  by  the  system,  (based  on  the  quad¬ 
rature  addition  of  terms  and  obtained  via  program  in  .Appendix  E,  M) 
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detected  range  is  limited  by  the  system  (quadrature  addition) 


90 

1.4e+004 


height  of  target  above  surface*3m  height  of  sensor  above  surface=3m 


Figure  Zjb  The  maximum  deteaed  range  is  limited  to  7,000m  due  to  the  method,  instead 
the  12.j65m  which  is  the  geometric  horizon  (based  on  the  quadrature  addi¬ 
tion  of  terms  and  obtained  via  program  in  Appendix  E.M) 
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VI.  GENERALIZATION  OF  THE  RESULTS  FOR  FIXED  SENSITIVITY 


A.  MOTIVATION 

As  seen  in  the  previous  chapter  the  condition  for  fixed  sensitivity  leads  to  polar 
curves  for  maximum  triangulation  range  using  both  a  LBL  and  a  SBL.  The  crossing  point 
for  these  curves,  Ocross’  defines  the  complimentary  application  domains  for  the  respective 
baselines.  In  this  chapter  some  simple  guidelines  for  predicting  Bcross-  maximum  SBL 
range,  (at  6=0°),  and  the  maximum  LBL  range,  (at  6=90°),  are  derived.  This  is 
done  for  the  quadrature  addition  of  terms  and  the  more  conservative  sum  of  absolutes 

B.  A  SIMPLE  FORMULA  FOR  THE  LBL  &  SBL  CROSSOVER 
1.  Analysis  for  the  sum  of  the  absolutes 

It  is  observed  that  in  the  crossover  points  of  the  LBL  and  SBL  polar  plots  the 
normalized  range  in  respect  to  LBL  for  both  baselines  has  to  be  equal.  That  means 


Hul  ^  Rgi 
LBL  "  LBL 


(70a) 


Taking  into  account  (61),  and  (62),  (70a)  becomes 


(70b) 


69 


or 


sin‘(0) 


=  COS' 


,|c-^)'+(l+fl^V 


(1+0*)- 


sm. 

LBL 


(70c) 


Furthermore  (70c)  becomes 


tan-(0) 


(71) 


Normally  in  (71),  the  terms  of  the  form  )‘  are  much  bigger  than  the  terms  of  the  form 
( 1 -t-  ''  )■  and  dominate,  so  (71)  can  be  written 

For  the  typical  case  when  c^=c^=t'.  Equation  (72)  becomes 

tan(e)=(^]  (73a,  b) 


and  finally 


Across  =  Otan  |  ^ 


(74) 


It  is  worth  noting  for  the  correction  of  (70)  that  if  and  c^=c^  then  (74)  is  exact 
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2.  Analysis  based  on  the  quadrature  addition  of  terms 


Similarly  to  the  previous  case  at  the  crossover  points  the  normalized  ranges  of 
both  baseline  systems  are  equal.  This  implies  that 

=  (75) 


and  finally  after  (63)  and  (65),  (75)  gives 


(c^)  Mn^(e'-)-0  5cos7(9^) 


+  sin 


•  .  f  (<r')^$in^(8*)-0  5cos7(9*)  ,  ■>  -o  W  cb; 

in‘(0^)  =  +  sin-(e^) 


(76) 


and  after  (24),  and  (25),  (76)  become 


|f^)‘sin‘(fl)-0  5cos7(9)  .  2/n\  f  (c^)*cos*(9)-0  5sm^(9) 

- 1 - +  sm  (0)  =  l - ^ - +  COS 


(77) 


or 


(c^^m^9W)5(l-5m=(9)]  ,  f  5sm7(9) 

- ; - +  sin‘(0)  = 


^1  -sin‘(0) 


(Sr  I  (78) 


After  assuming  that  c‘^=c^=c  and  substituting  sin^(e)  with  "y,"  in  (78),  it  follows  that 


0.5c*>'-0.25+  1.25_>'  =  (0.5c' -0.5c‘>'+  1-1. 


(79a) 
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and  finally 


0.5c-  +  1.25  +  0.5c-f^ 


+  1.25 


5  +  0,5r|  = 


(7%) 


Solving  for  "y,"  that  is  sin'(0),  and  after  back  substitution  into  (77),  it  follows  that 


sin"(6) 


0 

l2MOSc=-K)5c=(^)'+125(§f)' 


(80) 


It  can  be  observed  at  this  point,  that  after  dividing  both  terms  of  the  ratio  by  c'  and  ne¬ 
glecting  the  very  small  terms  that  follow,  (78)  becomes 


sin‘(0) 


(81a) 


or 


fai) 

sin  (6)^  (81b) 


and  finally 


r 


Across  =osin 


=  atan 


SBl. 

LBL 


(82) 
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3.  Observations 


After  the  analysis  of  both  methods  have  been  completed  the  very  obvious  result 
is  that  the  crossover  point  depends  only  on  the  ratio  of  the  two  baselines  and  it  is  approxi¬ 
mately  independent  of  the  range  sensitivity  factor  c.  To  demonstrate  the  above  statement 
the  more  exact  (80)  and  the  approximate  (82)  crossover  angle  is  plotted  in  Figure  24  for  a 
fixed  baselines  ratio  SBL/LBL=0.2  versus  the  variable  normalized  range  sensitivity  The 
approximate  form  (82)  is  nearly  equivalent  to  (80). 

C.  SIMPLE  SCALING  FORMULA  FOR  DUAL  BASELINE  MODEL  RANGES 
1.  Analysis  for  the  sum  of  the  absolutes 

As  seen  in  Chapter  IV  the  sensitivity-fixed  triangulated  ranges,  normalized  by 
the  LBL,  become  maximum  at  the  angle  of  90“,  measured  with  respect  to  each  baseline. 
This  maximum  normalized  range  can  be  evaluated  from  expressions  (61)  and  (62)  for  the 
long  and  short  baselines  respectively,  as  follow 


-JJJ— l0t=9(r  - 


(1+a^)' 


(83a) 


and 


I  _  I  f  SBL  y 


After  assuming  that  c^=c^=c,  and  f/=a^=l,  (81)  becomes 


(83b) 
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LHL  “ 


and 


(84a,  b) 


2.  Analysis  for  the  quadrature  addition  of  terms 

Following  the  same  procedure  described  above  (63),  and  (64)  give 


^ lAlunn)  I 
LBL 


(85a) 


and 


Finally  supposing  c‘-=cf=c,  (83)  follows  that 


(85b) 


(86a,  b) 


As  expected,  a  factor  of  2''^  increase  in  the  less  conservative  quadrature  addition  of  the 
terms  is  seen  by  comparing  with  (84). 

3.  Observations 

The  most  obvious  result  of  the  analysis  above  is  that  the  scale  of  the 
curves  increases  linearly  with  the  normalized  range  sensitivity.  More  specifically  for  the 
sum  of  absolutes  paradigm  the  maximum  ranges  satisfy  the  following  rules  of  thumb 
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=  T  xLfl/4m)  and  RsB!(m)  =  ^'x.SBL(m) 


(87a,  b) 


Similarly,  for  the  quadrature  addition  of  terms  paradigm  the  maximum  ranges  satisfy  the 
following  relations 


Ruii.{m)  =  -^xLBLim)  and  RsHL{m)  =  -^xSBL{m)  (88a,  b) 


For  both  paradigms  the  crossover  angle  is  given  by  the  rule  of  thumb 


Across  =  <?tan 


(89) 


The  general  characteristics  of  the  dual  baseline  triangulated  ranges  and  baseline  crossing 
points  are  shown  on  Figure  25. 


75 


angle  theta  In  degrees 
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comparison  between  crossover  angle  approximations 
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Figure  24b:  Crossover  angle  of  long  and  shon  baseline  based  on  the  accurate  expression 
(80)  and  the  approximate  expression  (82)  The  discrepancy  between  the  two 
expressions  becomes  obvious  as  the  scale  on  the  y  a.xis  goes  to  the  third  deci 
mal  point,  (via  program  in  .Appendix  E.  N) 


case  A 


(case  A  sum  of  absolutes  ,  case  B  .  quadrature  addition  of  terms  ) 


Figure  25  :  Properties  which  characterize  the  performance  of  the  dual  baseline  triangula¬ 
tion  scheme  The  maximum  normalized  range  which  can  be  detected  by  the 
system  is  linearly  dependent  on  the  sensitivity.  The  crossover  angle  of  the 
long  and  short  baseline  approximately  depends  only  on  the  ratio  of  the  two 
baselines 
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VII.  CONCLUSIONS 


A.  RESULTS  OF  THE  WORK 

A  dual  baseline  scheme  based  on  the  triangulation  principle  proposed  for  passive 
range  measurements  was  presented  and  its  performance  was  examined  in  this  work  In 
Chapter  II  the  approximations  on  which  the  proposed  model  was  based  were  examined, 
and  the  general  characteristics  of  the  triangulation  principle  were  reviewed  The  conclu¬ 
sion  was  that  the  approximations  R»BL,  and  R=R,  =R2  have  not  noticeably  changed  the 
accuracy  of  the  model,  and  that  for  angles  close  to  0°  and  1 80“  due  to  "geometric  dilution" 
no  range  measurements  can  be  performed.  Furthermore  it  was  shown  that  the  precision  in 
the  range  measurement  was  dependent  on  the  precision  in  the  bearing  measurement.  Lastly 
it  was  shown  that  the  more  distant  the  target  to  be  triangulated  is  from  the  baseline  the 
more  inaccurate  the  measurement.  In  Chapter  III,  feasible  modeling  for  a  dual  baseline  tri¬ 
angulation  scheme  was  examined.  It  was  shown  that  two  orthogonal  baselines  eliminate 
the  "geometric  dilution"  problem  and  that,  in  general,  the  larger  baseline  will  have  a 
higher  measurement  accuracy.  Finally,  for  further  study  of  the  characteristics  of  the  dual 
baseline  scheme  two  gauges  were  proposed.  The  first  one  is  based  on  the  quadrature  addi¬ 
tion  of  the  terms  which  express  the  range  sensitivity  of  each  sensor  for  each  baseline,  and 
the  second  one  is  approximate  and  is  based  on  the  sum  of  absolutes  of  the  same  terms  It 
was  shown  in  Appendix  D  that  the  latest  one  is  more  conservative  with  respect  to  the  er¬ 
ror  in  the  range  measurement,  and  it  can  be  used  for  arbitrary  bearing  precision  of  the 
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sensors.  Instead,  the  less  conservative  quadrature  addition  of  terms  can  be  used  only  for 
the  case  when  the  sensors  have  the  same  bearing  precision  In  Chapter  IV  the  perform¬ 
ance  of  the  dual  baseline  model  was  studied  by  the  above  mentioned  gauges  through  the 
polar  plots  of  the  normalized  range  versus  the  target  orientation  angle.  The  conclusion 
from  this  chapter  was  that  with  the  use  of  the  conservative  sum  of  absolutes  the  predicted 
maximum  range  is  2' '  less  than  the  range  predicted  by  the  quadrature  addition  of  terms  In 
Chapter  V  the  physical  limitation  in  the  range  measurements  due  to  optical  horizon  was 
considered.  The  effect  of  refraction  introduces  a  complicated  situation  in  the  prediction  of 
the  optical  horizon.  Thus,  a  more  conservative  estimation  for  the  optical  horizon  based  on 
the  geometrical  consideration  alone  was  taken.  It  was  shown  in  this  chapter  that  the  pre¬ 
dicted  maximum  range  depends  more  generally  on  sensor  and  target  height  above  the 
earth's  surface,  and  the  ratio  of  the  two  baseline  lengths.  Therefore  the  range  measurement 
can  be  limited  either  due  to  the  optical  horizon  or  the  dual  baseline  system  In  Chapter  VI 
the  general  features  of  the  polar  plots  that  are  characteristics  of  the  dual  baseline  method 
were  considered.  It  was  shown  that  the  maximum  range  scales  with  the  range  sensitivity  of 
the  baselines.  In  the  case  when  the  normalized  sensitivity  in  both  baselines  is  the  same 
then  the  general  appearance  of  the  two  curves  will  depend  only  on  the  ratio  of  the  long 
baseline  to  the  short  baseline  More  specifically,  the  crossover  angle  of  the  two  curves,  to 
a  high  degree  of  approximation,  depends  only  on  the  ratio  of  the  two  baselines,  and  it  is 
independent  of  the  range  sensitivities. 
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B.  FUTURE  DEVELOPMENTS 


The  dual  baseline  model  was  studied  under  the  assumption  that  the  mean  bearing  er¬ 
ror  was  zero,  but  in  the  real  world  that  is  not  true.  Thus  future  work  is  needed  to  enhance 
the  model  to  handle  systematic  errors  in  the  bearing  measurements.  This  future  work 
could  employ  the  "circulation  method"  instead  of  triangulation.  A  second  point  for  future 
development  could  be  the  implementation  of  the  dual  baseline  model  with  some  "smart 
electronics"  which  would  switch  between  baselines  at  a  performance  crossover  point  , 
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APPENDIX  A:  LIST  OF  SYMBOLS 


A.  SYMBOLS  IN  ENGLISH 

BL:  the  baseline  (platform  length)  in  the  triangulation  method  (m) 
c.  normalized  range  sensitivity  in  triangulation  method 
normalized  range  sensitivity  in  long/short  baseline 
d:  distance  to  the  geometrical  horizon  (m) 

LBL:  the  long  baseline  in  the  dual  baseline  triangulation  scheme  (m) 

R;  the  range  from  the  middle  of  the  baseline  to  the  target  (m) 

R^.  radious  of  the  earth  (6.371  x  10®  meters) 

Rlpl  the  range  from  the  middle  of  the  long  baseline  to  the  target  (m) 

Rjgj^:  the  range  from  the  middle  of  the  short  baseline  to  the  target  (m) 

Ri,,:  the  range  from  the  first  second  sensor  to  the  target  (m) 

SBL:  the  short  baseline  in  the  dual  baseline  triangulation  scheme  (m) 
s:  height  of  the  sensor  above  the  surface  of  the  earth  (m) 
t:  height  of  the  target  above  the  surface  of  the  earth  (m) 

B.  SYMBOLS  IN  GREEK 

a  :  ratio  of  the  mean  values  of  the  bearing  errors  in  the  triangulation  method 
:  ratio  of  the  mean  values  of  the  bearing  errors  in  long/short  baseline 
A6„;  :  bearing  error  of  the  first/second  sensor  (rad) 
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6  :  the  target  orientation  angle  measured  from  the  middle  of  the  baseline  (degrees) 

6  ‘  ^  angle  to  the  target  measured  from  the  middle  of  both  baseline  (degrees) 

6i ;  angle  to  the  target  measured  from  the  first/second  sensor  (degrees) 

X|  independent  random  variable  having  zero  mean  and  unity  standard  deviation  in 
long/short  baseline  in  respect  of  first/second  sensor 
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APPENDIX  B:  MATHEMATICAL  PROOFS 


A.  SINGLE  BASELINE  MODEL 
1.  Single  baseline  model  analysis 

The  distance  from  the  middle  of  the  baseline  to  the  target  given  by  Equation 
(4)  is  repeated  here 


R  =  ^ 


sindi 


1-  frsme: 


(B.  1) 


The  derivatives  of  (B. ) )  with  respect  0,,  and  0,  are  respectively 


■ 

/ 

“ 

cosHf 

1 

l-(fsine:]  - 

sinBismBj 

-^(01.0;./?) 

sin(0;-0i( 

v 

1 

(B.  2) 


dF  _ 


( ]  sinBisinfl; 


+g(0,,B2.R) 


(B  3) 


where 


g(8i.e:,R)  = 


( 


sin*(0:-0| ) 


sinB; 


+  sin0 


(B  4) 
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and 


(B.5) 


The  expressions  (B.2),  (B.3),  and  (B.  4)  above  are  approximate  because  the  derivatives  of 
the  small  terms  proportional  to  (BL/R)'  on  the  RHS  of  Eq.  (B.  1)  have  not  been  included 
The  terms 


Sin  (02-6) ) 


and 


sin  (02”^  I ) 


(B  6) 


can  be  factored  out  of  Equations  (B  .2),  and  (B  4).  Also  the  terms 


sin(02-0|  J 


and 


sin(92-Bi) 


(B  7) 


can  be  factored  out  between  Equations  (B  3),  and  (B.4),  Then  Equation  (B  2),  after  Equa¬ 
tion  (B  6),  becomes 


(B.  8a) 


where 


fi  = — — (cosei  +sin(Jicot(e;-ei)) 


(B  8b) 
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and 


SinB; 


h 


sin(H2-ft|) 


cot  (02  -0l)- 


sm  Hi  COS  Hi 


(B.  8c) 


Similarly,  Equation  (B.3),  after  Equation  (B  7),  becomes 


dR 


(B.  9a) 


where 


)~  ~ sin  (9:  -  e I )) 


(B  9b) 


and 


sin 


001(02 -0|)  + 


V 


( Is  ) 

l-(fismB:] 


J 


(B,  9c) 


Equation  (B.  5),  after  (B.  8),  (B  9),  and  after  dividing  by  (B.  1),  becomes 
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^sinB:]  (cosfli4«m9|Cot(93-Bi)>+sin03 
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(cos0:-sinB2Coi(B:— 0  ))-sin0 


COt(02— 0j  )-4— 


Sind  I 


■Minfli  Jl-(^s 


(Adz)  (B  10) 


as  cited  in  the  main  text. 

2.  Single  baseline  model  analysis  (R»BL) 

A  reasonable  assumption  that  R»BL  can  be  made,  and  it  follows  from  Equa¬ 
tion  (4) that 


R 


BL  sinSj+sinSi 
2  sin(02-di) 


) 


The  derivatives  of  (B.  1 1 )  with  respect  0,,  and  6,  are 


BL  I  cos9i 
2  ,  sin  ) 


^(61.62)) 


and 


dR 


BlI  cos  9; 

2  I  sm(9j-e,) 


(B.  11) 


(B  12) 


(B.  13) 
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where 


^(e,.0;) 


sin"(B:-8j) 


(B.  14) 


Equations  (B.  12),  and  (B.  13)  after  (B.  14)  can  be  re-expressed  as  follow 

^  =  +C0t(e2 -BiXsinBi  +sine:))  (B.  15) 


and 


dR 

cm2 


-  0 1  )(sin  0:  -  sin  0 : )) 


(B  16) 


Combining  Equations  (B.15),  and  (B.16),  and  after  dividing  by  (B.l  1)  leads  to 

f  =  (dS:5:+“'<9:-e'))  w9.)+(=fS5r-cot(e:-e,))  (4e,)  (b,  h) 

as  cited  in  the  main  text. 

B,  DUAL  BASELINE  MODEL 

Following  the  substitutions  suggested  by  Equations  (23),  (24),  (26),  and  (27)  the 
expressions  for  obtained  for  the  various  approximation  cases  of 

the  dual  baseline  model  analysis  introduced  in  Chapter  III. 
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1. 


Dual  baseline  model  analysis 


a.  Long  baseline  (LBL) 
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f.-(^.u,8‘)  = 

(cosfl;+sin0icol(8;-e{  ]  j-sinef^l-f^sinfl;  )  coi^Hi-Hf 


+5m0 


(B  18) 


X:« 


6.  baseline  (SBL) 


^  J_ 
AH 


^  1  ^  sin ef  j ‘  [cos 0f+sin ef cot (Sj-ef  ]  )+sin efji- 

/  ...  nSV 

/  \  *  i  *10  07  to*  S'? 

...f^S  „S)  2JI  1  ' 

1  2/t  j 

smdi^  I-^^sinSf  j  45in6 

xf+ 


I 


J  )  (cos0f+5inflfcol| 

8f-0f))-sm8f^l-(^su.0f)' 

CO 

„s) 

sinBj^ 

r 

1- 

(^sm8f)’‘H,u,8fJl-(^smef) 

(B.  19) 


2.  Dual  baseline  model  (R»BL) 
a.  Long  baseline  (LBL) 


( V = r 

\  ff  AH  j  l, 


+cot  (02 -et 


cos  0* 


■cot 


(B  20) 
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Short  baseline  (SBL) 


f  cos 


V  sin  Hx-«in  H‘| 


+  col  0 


xf+ 


sinBA+sinH', 


-cot  0-  -0; 


(B.2!) 
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APPENDIX  C:  EXACT  EXPRESSION  FOR  THE  RANGE 

A.  GENERALITIES 

From  Figure  2  the  expression  (3b),  which  relates  the  range  of  the  target  to  the  apex 
angles  on  the  baseline,  had  been  derived.  Solving  (3b)  for  R  gives  the  exact  expression 
that  has  to  be  used  to  solve  the  problem  in  the  triangulation  method. 

B.  ANALYSIS 

In  order  to  simplify  the  notation  the  following  definitions  are  introduced 

ATI- H  sin  (02  -01 ) 

^  s  sin  02 
sw  =  sin0i 

Then  from  (3b) 


After  substituting 

Z=(f)  (C.3) 


(C.  la) 
(C  lb) 
(C.  Ic) 
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symbolic  processing  using  MAPLE  [Ref.  1 9],  produces 


4Z‘^.si  ‘(-l  +Z^sy^)  =  (Z~sy^  -IZ'^.sy^sM’^  +Z*.sit  ‘  -.sn  ')'  (C.  4) 


and  then  again  symbolic  processing  produces  the  two  solutions 


2  j  -Isx^sy^  - Isw^sx^  +  4  +sx^) 

2  -As\i’^ sy^ sx^  -  Ay'*  +  IsM'^sy^  -  5^“* 


(C  5a) 


-  1  -2sX^5V^  -  25M'^SX^  -  4  J-SX^Sy^SH'^(-  \  +5X^) 

(Z2)2  =  i- - ^ - 


2  -Asw^sy^sx^  —  +  2^^^.^^  -  w'* 


(C.  5b) 


It  can  be  shown  using  MAXSYMA  [Ref  20],  that  the  root  (C.  5a)  is  real  and  (C  5b)  is 
imaginary.  Therefore  the  exact  expression  for  normalized  range  in  terms  of  angles  0,.  and 
0,  is 


— 45M'^^^sx^  —sy*  +  Isw'^sy'^  -  5M’'’ 
-Isx^sy^  —  Isw^sx^  +  4  J-sx'^sy^sM'^i-l  +sx^) 


(C.6) 
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APPENDIX  D:  RANGING  ERROR  IN  TRIANGULATION  METHOD 


A.  GENERALITIES 

In  single  triangulation  the  target  lies  at  the  intersection  of  the  two  lines  that  run  from 
the  sensors.  These  two  lines  make  angles  with  the  baseline  that  vary  due  to  mechanical 
limitations  of  the  sensors.  The  variation  in  the  angles,  A0,  and  A6,.  determine  the  varia- 
*'on  in  the  measured  range  of  the  target. 

B.  APPROXIMATE  RANGING  ERROR  ANALYSIS 

Figure  26  shows  the  geometry  of  range  measurement  in  the  single  baseline  triangula¬ 
tion  scheme.  The  quadrangle  KLMN  is  established  due  to  variation  in  the  angles.  The  lines 
DE,  and  FG  have  values  AR,  and  AIL  respectively,  and  approximately  represent  the  un¬ 
certainties  in  the  range  of  the  target  from  the  sensors.  The  line  MK  represents  the  total  un¬ 
certainty  due  to  the  system  and  it  has  value  AR.  In  first  approximation  the  error 
quadrangle  can  be  treated  as  a  parallelogram  and  the  following  relation  is  derived  from 
parallelogram  on  Figure  26 


ARq  =  ARparauelogram  =  ARp 


(D  1) 


or 


A/?,.  =  ^(Ay?,)-  +  (A/?2)‘  +  2|A^,|  •  |A/?;|cos(t) 


(D  2) 
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where  "(f)"  is  the  angle  between  AR,  with  AR^. 

For  R»BL,  as  shown  in  Figure  27,  <()  =  0‘,  and  therefore 

cos(t)  =  1  (D  3) 

Thus,  it  follows  from  Equation  (D.  2)  that 

Mp-M<lM,l  +  lM2l  (D.4) 

and  therefore  via  comparison  with  (D.2)  it  follows  that 

<fA/t,i  +  lA/i:l  (D.5) 

In  fact,  under  the  case  that 


|AR,|  =  |AR:| 


(D.6) 


e  g  at  a  point  of  symmetry,  and 


COS({)=  1 


(D7) 


the  sum  of  quadratures  estimate  will  be  2''^  smaller  than  the  sum  of  absolutes.  So,  in  gen¬ 
eral,  the  sum  of  absolutes  is  a  more  conservative  gauge  for  the  lack  of  precision  in  the  tri¬ 
angulation  measurment. 
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ERROR 


QUADRANGLE 


BL 


Figure  26:  Geometry  for  ranging  error  in  single  baseline  triangulation 
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APPENDIX  E:  MATLAB  CODES 


A.  SINGLE  BASELINE  MODEL 

clear; 

b=input  ('put  the  ratio  RyBL='); 

y=l./b; 

a=input('put  the  ratio  DTH2./DTH1-); 
step=pi./40; 
theta=0:step.pi; 
u=sin(theta); 
v=cos(theta); 
q=tan(theta); 

j=i  /q; 

c=(b.^2+.5/^2+b.*v)/'.5;  %  c=R2/BL  eq.  9a  ; 

d=(b.^2+.5/'2-b.*v)/'.5;  %  d=Rl/BL  eq.  9b  ; 

thetal=asin(u  *b./c); 
e=sin(thetal); 
f^cos(thetal); 
theta2=asin(u.*b./d); 
g=sin(theta2); 
h=cos(theta2); 
theta21=asin(l./d.*e); 
z=tan(theta2 1 ); 
k=l./z; 
l=(,5.*y), 
m=l.*g; 

n=(I-ni.^2).^.5; 
p=l,*e; 

q=(l-p.^2).-5; 
r=g.*q+e,*n; 
w=k-(1.^2.'»e,*f./q.^2); 
x=k+(l,''2.*g.*h./m.^2); 

B=n,  *(f+e.  *k)+g.  *q.  '^w; 

C=q  *(h-g.*k)-e.'*n.*x; 

D=(f/(g+e)+k); 

E=(h/(g+e)-k); 

F=(sqrt(l-(y.*u).^2)./(y,*u)), 

%  s  represents  the  general  solution  for  single  baseline  model  eq  7 
s=((B./r).^2+((C./r),*a).^2).^.5; 


% 

eq.  10b ; 

% 

eq.  10a ; 

% 

eq  11  ; 

%  1=BL/2R; 

%  m=(BL/2R)*sin(theta2); 
%  p=(BL/2R)*sin(theta  1 ); 
%  r=denom.  of  eq  .  7  ; 
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%  t  represents  the  approximate  solution  (R»BL)  eq  15  ; 
t=(D.^2+(E.*a)/'2).^.5; 

%  z  represents  the  approximate  solution  (R»BL,R1=R2=R)  eq  21  , 

z=((.5  *j+F)/'2+((.5.*j-F).*a).^2i/'.5; 

theta=theta.*180./pi; 

plot(theta,s,'*',theta,t,'o',theta,z,'+'); 

axis([0  180  0  . ]); 

xlabel(Theta  in  degrees'); 
ylabelCjdelta  R  /  (R  *  delta  theta)|'); 

titleCcomparison  of  single  baseline  model  sensitivity  curves'); 

gtext("^  for  more  accurate  case  (equation  7)'); 

gtext('o  for  R»BL  approximation  case  (equation  15)'), 

gtextC  +  for  R»BL  &  R=R1=R2  approximation  case  (equation  21)'); 

gtext('R/BL='); 

gtext('a='); 


B.  DUAL  BASELINE  MODEL 


clear; 

rovlbl=input('put  the  ratio  R/LBL-);  % 

y=l./rovlbl;  % 

rovsbl=input('put  the  ratio  R/SBL- );  % 

x=l./rovsbi,  % 


a=input('put  the  ratio  DTHL2/DTHL1-); 

aa=input('put  the  ratio  DTHS2./DTHS1=); 

step=pi./40; 

theta=0;step.pi; 

u=sin(theta); 

v=cos(theta); 

m=(  1  -sin(theta)/'2)/'.  5 , 


rovlbl=R/LBL; 

y=LBL/R; 

rovsbl=R/SBL; 

x=SBL/R; 


c=y  *u; 


b=(l-c/'2)/'.5, 

d=(l-(-v.^2))/'‘;; 
g=x  "^(-v); 
Hl-g/'2).^.5; 


k=v./u,  %  k=cot(theta); 

%  c=(LBL/R)'''sin(theta); 

%  b=(  1-((LBL/R)'*sin(theta)''2)^.5; 

%  d=(J-(-cos(theta)^2)'',5, 

%  g=(SBL/R)*(-cos(theta)); 

%  Hl-((SBL/R)*(-cos(theta)))''2)^.5; 
%  the  following  relation  is  for  LBL  model  eq  33a  ; 
l=((l/2  *k+b./c)/'2+((l/2.*k-b./c)  *a)/'2)/' 5, 

%  the  following  relation  is  for  SBL  model  eq  33b  ; 

s=((l/2.*d/(-v)+f/g)/'2+((l/2  *d./(-v)-f./g),*aa)/'2)/'.5; 

theta=theta.  *  1 80./pi; 

plot(theta,l,'o',theta,s,"*'); 

axis([0.  180.  0,  2.5*10^4]); 

titleC  dual  baseline  model  sensitivity  curves'); 
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xlabeiCTheta  in  degrees'); 
ylabelCDelta  R/(R  *  Delta  Theta'); 
gtextCo  for  long  baseline'), 
gtextC*  for  short  baseline'), 
gtextCRTLBL- ), 
gtext('R/SBL- ); 
gtext('a='); 

gtext('ratio  of  bearing  errors  in  long  baseline='), 
gtext('ratio  of  bearing  errors  in  short  baseline^'), 

C.  RANGE  LIMITATIONS  FOR  FIXED  TARGET  ORIENTATION 

clear; 

maxrovbl=input('  input  the  maximum  ratio  of  r/bl  '); 
minrovbl=input('  input  the  minimum  ratio  of  r/bl  '), 
thetv=l:3; 

thetv(  1  )=  input  ('  angle  1  in  degrees  other  than  0, pi'); 
thetv(  1  )=pi/ 1 80*thetv(  1 ), 

thetv(2)=  input  ('  angle2  in  degrees  other  than  0,pi'); 
thetv(2)=pi/l  80*thetv(2); 

thetv(3)=  input  ('  angle3  in  degrees  other  than  0,pi'), 

thetv(3)=pi/180*thetv(3); 

nmax=input('  input  the  #  points  '); 

long=l  nmax; 

longl  =  l  :nmax; 

long2=l  :nmax; 

long3=l  :nmax, 

longv=long'*thetv; 

difr=  (maxrovbl-minrovbl)/(nmax-l); 

for  k=l:3, 

thet=thetv(k); 

for  i=l:nmax; 

rovbl(i)=  minrovbl+difr*(i- 1 ); 
rovbll=rovbl(i); 
v=sin(thet) ; 
z=tan(thet); 

x=l./z; 

long(i)=  l/2.*x  +  (sqrt(l.((l  /rovbll).*v)/'2))./((l./rovbll)  *v)) ;  %eq  47 

longv(i,k)=long(i); 

end  ; 

end; 

for  i=l  :nmax  ; 
longl(i)=longv(i,l); 
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long2(i>="!.';ngv(i.2). 

lonp’'^-;-longv(i,3), 

fi' 

plot{rovbl,long  1  ,'o',rovbl,long2,'*’,rovbl,long3.'x')  . 

title  ( '  symmetrical  analysis  ' ) 

xlabel  ('  range/base-line  ')  ; 

ylabel  ('  delta  R/(R*  delta  theta)  ') , 

gtextCo  for  target  approching  at  an  angle  --  degrees’); 

gtextC*  for  target  approching  at  an  angle  —  degrees'); 

gtextfx  for  target  approching  at  an  angle  --  degrees'); 


D.  LONG  BASELINE  POLAR  PLOT  (SUM  OF  ABSOLUTES) 

clear, 

%al=input('put  the  ratio  of  the  bearing  errors  DTH2L/DTH 1  L='), 
al=l; 

%cl=input('put  the  critical  value  of  the  normalized  sensitivity^'); 
cl=500, 

theta=0  0  05:2*pi, 

%  rlovlbl  represents  the  eq.55d  , 

rlovlbl=abs((sin(theta).^2  '^((cl/'2+(l+al)/'2.)/(l+al)/'2))/'  5); 
polar(theta,  rlovlbl), 

titleCnormalized  range  detected  by  long  baseline  (sum  of  absolutes)'), 
gtextCratio  of  the  bearing  errors  in  long  baseline='), 
gtext('normalized  sensitivityof  the  long  baseline=’); 
gtext('detected  range'); 

E.  SHORT  BASELINE  POLAR  PLOT  (SUM  OF  ABSOLUTES  ) 

clear; 

%as=input('put  the  ratio  of  ratio  of  the  bearing  errors  DTH2S/DTH1  S='); 
as=l, 

%f=input('ratio  of  the  two  baselines  SBL(m)/LBL(m)='), 

1^0.2, 

%  cs=input('put  the  critical  value  of  the  normalized  sensitivity^'); 
cs=500, 

theta=0:0,05;2*pi, 

%  rsovlbl  represents  the  eq  55f ; 

rsovlbl=abs((cos(theta)/'2  *((cs/'2+(l+3s)/'2)  /(l+as)/'2) 
polar(theta,rsovlbl); 

titleCnormalized  range  detected  by  the  short  baseline  (sum  of  absolutes)'); 
gtext(’LBL=5  x  SBL'), 

gtextfratio  of  the  bearing  errors  in  short  baseline='); 
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gtext('normalized  sensitivity  of  the  short  baseline-), 
gtext('detected  range'); 


F.  LBL  &  SBL  CARTESIAN  PLOT  (SUM  OF  ABSOLUTES) 

clear, 

%cl=input('critical  value  of  the  normalized  sensitivity  in  long  baseline='), 
cl=500, 

%al=input('ratio  of  the  bearing  errors  in  LBL  (DTH2/DTH1)='); 
al=l; 

%f==input('ratio  (SBL(m)  /LBL(m)='); 
f^0,2; 

%cs=input('critical  value  of  the  normalized  sensitivity  in  short  baseline  ='); 
cs=500; 

%as=input('ratio  of  the  bearing  errors  in  SBL  (DTH2/DTH1  )='), 
as=l; 

theta=0:0.05:2*pi; 

rlovlbl=abs((sin(theta)/'2.*((cl.^2+(l+al).^2)./(l+al)/'2))/'.5);  %  eq  55d  ; 

rsovlbl=abs((cos(theta)/'2  *((cs.  ^2+(  1 +as)/'2)./(  I +as) /'2).  * . . 
f''2).^.5);  %  eq  55f , 

rmax=max(rlovlbl,rsovlbl); 

%polar(theta,rmax), 
theta=theta  *  180  /pi; 

plot(theta,rlovlbl,'.',theta,rsovlbl,'-.',theta,rmax,'o'); 

%gtext('al=2'); 

%gtext('cl=500'); 

%gtext('SBL/LBL=0.2'); 
xlabeK'theta  in  degrees'); 

%ylabel('rlbl/lbl,  rsbl/lbl,  rmax/lbl'); 

title(’cartesian  plot  of  the  range  detected  by  the  dual  baseline  system  (sum  of 
absolutes)'); 

gtext('ratio  of  the  bearing  errors  in  both  baselines^  1 '); 
gtext('normalized  sensitivity  in  both  baselines=500'); 
gtext('LBL=5  x  SBL'), 
gtext('range  detected  by  short  baseline'); 
gtext('range  detected  by  long  baseline'); 
gtext('range  detected  by  the  system'); 

G.  LBL  &  SBL  POLAR  PLOT  (SUM  OF  ABSOLUTES) 
clear, 

%cl=input('critical  value  of  normalized  sensitivity  in  long  baseline='); 
cl=500. 
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%a!=input('ratio  of  the  bearing  errors  in  long  baseline  (DTH2L/DTH1L)='); 
al=l; 

%f=input('ratio  (SBL(m)  /LBL(ni)- ); 
f=0.2, 

%cs=input('critical  value  of  normalized  sensitivity  in  short  baseline  ='), 
cs=500, 

%as=input('ratio  of  the  bearing  errors  in  short  baseline  (DTH2S/DTH1S)='); 
as=l; 

theta=0;0.05:2*pi; 

rlovlbi=abs((sin(theta)/'2.*((cl.^2+(l+al)/'2)./(l+al)/'2))/'.5);  %  eq.  55d  ; 

rsovlbl=abs((cos(theta)/'2.  *((cs/"2+(  l+as)/'2)./(  1  +as)/'2).  . 

f^2).^.5);  %  eq.  55f ; 

nnax=max(rlovlbl,rsovlbl); 

polar( theta, rlovlbl, ' . theta, rsovlbl,'-  ',theta,rmax,'o' ) ; 

title('range  detected  by  the  dual  baseline  system  in  any  direction(sum  of  absolutes)'), 
gtext('ratio  of  bearing  errors  in  both  baselines=  1'); 
gtext('normalized  sensitivity  in  both  baselines=500'); 
gtext('LBL-5  x  SBL'); 

H.  LONG  BASELINE  POLAR  PLOT  (QUADRATURE  ADDITION  OF  TERMS) 

clear; 

theta=0:0.05:2*pi; 

%cll=input('critical  value  of  the  normalized  sensitivity  in  long  baseline='); 
cll=500; 

%  rlovlbl  represents  eq.  61a  ; 

rlovlbl=abs(((  1  /2.  *(cll.  ^2 ,  *sin(theta).^2-0. 5 .  *cos(theta).  ^2)+sin(theta).  ^2)).  5 ); 

polar(theta,  rlovlbl); 

titleCnormalized  range  detected  by  long  baseline  (quadrature  addition  of  terms)'); 
gtext('ratio  of  the  bearing  errors  in  long  baseline=r); 
gtext('normalized  sensitivity  of  the  long  baseline=500'); 
title('detected  range'); 

I.  SBL  POLAR  PLOT  (QUADRATURE  ADDITION  OF  TERMS) 
clear; 

theta=0:0.05:2’'‘pi; 

%css=input('critical  value  of  the  normalized  sensitivity  in  short  baseline='); 
css=500; 

%f=input('ratio  SBL(m)  /LBL(m)- ); 
f=0.2; 

%  rsovlbl  represents  eq  61b, 

rsovlbl=abs(((l/2.’*(css/'2.*cos(theta)/'2-0,5  *sin(theta)/'2)+ 
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cos(theta).'^2)/‘f/'2).^0.5); 

polar(theta,rsovlbl); 

gtext('ratio  of  the  bearing  errors  in  short  baseline=r); 
gtext('normalized  sensitivity  of  the  short  baseline=500'); 

title('normalized  range  detected  by  the  short  baseline  (quadrature  addition  of  terms'), 
gtext('LBL=5  x  SBL'); 
gtext('detected  range'); 

J.  LBL  &  SBL  CARTESIAN  PLOT  (QUADRATURE  ADDITION  OF  TERMS) 

clear; 

%cll=input('critical  value  of  the  normalized  sensitivity  in  long  baseline='); 
cll=500; 

%ff=input('ratio  SBL(m)  /  LBL(m)- ); 
fF=0.* 

%css=input('critical  value  ot  the  normalized  sensitivity  in  short  baseline='); 
css=500; 

theta=0;0.05:2*pi; 

%  rlovlbl  represents  eq.  61a  ; 
rlovlbl=abs(((l/2  *(cll/'2.*sin(theta)/'2-... 

0.5,*cosrtheta)/'2)+sin(theta).^2))/O.5); 

%rsovlbl  represents  eq.  61b; 

rsovlbl=abs(((l  /2.*(css.^2.*cos(theta).^2-0.5  .*sin(theta).^2)+, . . 

cos(theta).''2).*fF.^2).^0,5); 

rmax=max(rlovlbl,rsovlbl); 

plot(theta,rlovlbl,'.',theta,rsovlbl,'-.',theta,rmax,'o'); 
xlabel('theta  in  degrees'); 
ylabel('rlbl/lb!,  rsbl/lbl,  rmax/lbV); 

title('cartesian  plot  of  the  range  detected  by  the  dual  baseline  systerr.(quadrature 
addition)'); 

gtext('ratio  of  bearing  errors  in  both  baselines  =1'); 
gtext('normalized  sensitivity  in  both  baselines=500'); 
gtext('LBL=5  x  SBL'); 
gtextC  range  detected  by  short  baseline'); 
gtext('range  detected  by  long  baseline'); 

K.  LBL  &  SBL  POLAR  PLOT  (QUADRATURE  ADDITION  OF  TERMS) 
clear; 

%cll=input('critical  value  of  the  normalized  sensitivity  in  long  baseline='); 
cll=500; 

%fF=input('ratio  SBL(m)  /  LBL(m)- ), 
fF=0.2; 
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%css=input('critical  value  of  the  normalized  sensitivity  in  short  baseline='); 
css=500; 

theta=0:0.05:2*pi; 

%  rlovlbl  represents  eq.  61a  ; 

rlovlbl=abs(((  l/2.*(cll.^2.*sin(theta).^2-0.5  *cos(theta)/'2)+sin(theta)/'2))/^0.5). 

%  rsovlbl  represents  eq.  61b  ; 

rsovlbl=abs(((  1  /2 .  *(css,  ^^2 ,  *cos(theta).^2-0. 5 .  *sin(theta).  ^2 )+. . . 

cos(theta)/^2).*fr.^2).^0.5); 

rmax=max(rlovlbl,rsovlbl); 

polar(theta,rmax); 

theta=theta.  *  1 80. /pi; 

titleCrange  detected  by  the  dual  baseline  system  in  any  direction  (quadrature 
addition)'); 

gtext('ratio  of  bearing  errors  in  both  baselines=]'); 
gtext('normalized  sensitivity  in  both  baselines=500'); 
gtext('LBL=5  x  SBL'); 
gtext(’ass=r); 

gtext('range  detected  by  the  system'); 

LIMITATIONS  IN  THE  MEASURMENT  OF  THE  RANGE  (S.  A) 

clear; 

%  S.  A  is  for  the  sum  of  absoluets; 

%al=input('ratio  of  the  bearing  errors  in  the  long  baseline(  DTH2L/DTH1L)='); 
al=l; 

%cl=input('critical  value  of  the  normalized  sensitivity  in  long  baseline='); 
cl=500; 

%lbl=input('length  of  the  long  baseline(m)- ); 

lbl=  100/20; 

theta=0:0.05;2*pi; 

%  rl  represents  the  eq.  55d; 

rl=abs((sin(theta).^2  *lbl.''2.*((cl.^2+(l+al).''2)./... 

()+al).^2)).''0.5); 

%as=input('ratio  of  the  bearing  errors  in  short  baseline  (DTH2S/DTH 1 S)- ); 
as=l; 

%cs=input('critical  value  of  the  normalized  sensitivity  in  short  baseline='); 
cs=500; 

%sbl=input('length  of  SBL(m)='); 
sbl=20/4; 

%  rs  represents  the  eq.  55f; 

rs=abs((cos(theta).''2.  ’*sbl.''2.  *((cs.''2+(  1  +as).^2)./. . . 

(l+as)^2)).^.5), 

rmax=max(rl,rs); 
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%  re  represents  the  radius  of  the  earth  (6.371  x  10^6  meters); 
re=6.371*  10.^6; 

%s=input('height  of  the  sensor  above  surface  (m)-); 
s=10/3; 

%t=input('height  of  the  target  above  surface  (m)- ); 
t=10/3; 

%  d  represents  the  eq.  68; 
d=(2.*re.*t).'^.5+(2.*re.*s).'^.5; 

%  r  represents  the  geometrical  horizon; 
r=(d .  ^2 .  *cos(theta).  ^2+d .  ^2 .  *  sin(theta).  ^2).  5 ; 
if  d>rmax 
polar(theta,r); 
hold  on 

polar(theta,rmax); 
else 

polar(theta,rmax); 
hold  on 
polar(theta,r); 
end; 

title('detected  range  is  limited  by  the  geometric  horizon  /  by  the  system  (sum  of 
absolutes)'); 

gtext('long  baseline  length®  1 00/20m'); 
gtext('short  baseline  length  =20/4  m'); 
gtext('height  of  the  target  above  surface=I0/3m’); 
gtext('height  of  the  sensor  above  surface  =l0/3m'); 
gtextCgeometric  horizon'); 
gtext('range  detected  by  the  system'); 
gtext('detected  range'); 

M.  LIMITATIONS  IN  THE  MEASURMENT  OF  THE  RANGE  (Q.  A) 

clear; 

%  Q.A  is  for  the  quadrature  addition  of  terms; 

%cll=input('critical  value  of  the  normalized  sensitivity  in  long  baseline®'); 
cll=500; 

%lbl=input('length  of  the  long  baseline(m)='); 
lbl=  100/20; 

%sbl=input('length  of  the  short  baseline(m)- ); 
sbl=20/4; 

%css=input('critical  value  of  the  normalized  sensitivity  in  short  baseline®'); 
css=500; 

theta=0:0.05;2’'‘pi; 

%  rl  represents  eq.  61a; 
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rl=abs(((  1 /2.  *(cIl/'2,  *sin(theta)/‘2-0.5.  *cos(theta)/'2)+. . . 
sin(theta).^2).*lbl.^2).^0,5); 

%  rs  represents  eq.  61b; 

rs=abs(((l/2.*(css/'2.*cos(theta)/'2-0.5.*sin(theta)/'2)+... 

cos(theta).^2).*sbl.^2)/'0.5); 

rmax=max(rl,rs); 

%  re  represents  the  radious  of  the  earth  (6.37  lx  10^6  meters); 
re=6,371*  10,^6; 

%s=input('height  of  sensor  above  surface  (m)='); 
s=10/3; 

%t=input('height  of  terget  above  surface  (m)='); 
t=10/3; 

%  d  represents  eq.  68; 
d=(2.*re.*t).^.5+(2.*re.*s).^.5; 

%  r  represents  the  geometrical  horizon, 
r=(d .  ^^2 .  *  cos(  theta) .  '^2+d .  ^2 .  *  sin(theta).  ^2).  5 ; 
if  d>rmax 
polar(theta,r); 
hold  on 

polar(  theta,  rmax); 
else 

polar(theta,rmax); 
hold  on 
polar(theta,r); 
end; 

titleCdetected  range  is  limited  by  the  geometric  horizon/by  the  system  (quadrature 
addition)'); 

gtext('long  baseline  length=100/20m'); 
gtext('short  baseline  length=20/4m'); 
gtext('height  of  the  target  above  the  surface  =  10/3m'); 
gtext('heigt  of  the  sensor  above  the  surface=  10/3m'); 
gtextCgeometric  horizon'); 
gtext(’range  detected  by  the  system'); 
gtext('detected  range'); 

N.  CROSSOVER  ANGLE 
clear; 

%  f=input('ratio  of  the  two  baselines  SBL(m)/LBL(m)='); 
f=0.2; 

%  c  is  the  normalized  sensitivity  of  both  baselines  which  assumed  to  be  the  same; 

%  clower=  input  ('lower  limit  of  the  normalized  sensitivity='); 
clower=500; 
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%  cupper+input  ('upper  limit  of  the  normalized  sensitivity  ='); 
cupper=1000; 

%cstep=input  ('step  of  normalized  sensitivity'); 
cstep=  1 0; 

c=clower :  cstep  xupper; 

theta l=atan(f*(c./c));  %  eq.  82  ; 

theta2=asin(sqrt((0.25+0. 5 .  *c/^2 .  *f  "'2+. 

f^2),/(1.25+0.5.*c.'^2+0.5.*c.^2.*f^2+f^2)));  %  eq.  80  ; 

theta  1  =theta  1 .  *  1 80/pi; 

theta2=theta2.  *  1 80/pi; 

plot(c,thetal); 

hold  on; 

plot(c,theta2,'.'); 
axis([500  1000  11.3  11.4]); 
xlabel('normalized  range  sensitivity'); 
ylabel('angle  theta  in  degrees'); 
title('comparison  crossover  angle  approximations'); 
gtextCaccurate  cross  angle  (eq.  80)'); 
gtext('approximate  cross  angle  (eq.  82)'); 
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